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Abstract 


This  paper  i^esents  an  algorithm  for  the  generation  of  mechanical  assembly  sequences  and  a  proof  of  its  conectness 
and  completeness.  The  algorithm  employs  a  relational  model  of  assemblies.  In  addition  to  the  geometry  of  the 
assembly,  this  model  includes  a  representation  of  the  attachments  that  bind  one  part  to  another.  The  problem  of 
generating  the  assembly  sequences  is  transformed  into  the  problem  of  generating  disassembly  sequences  in  which 
the  disassembly  tasks  are  the  inverse  of  feasible  assonbly  tasks.  This  transformation  leads  to  a  decomposition 
tqjproach  in  which  the  problem  of  disassembling  one  assmbly  is  decomposed  into  distinct  subproblems,  each  being 
to  disassemble  one  subassembly.  It  is  assumed  that  exactly  two  parts  or  subassemblies  are  joined  at  each  time,  and 
that  whenever  parts  are  joined  forming  a  subassembly  all  cont»:ts  between  the  parts  in  that  subassembly  are 
established.  The  algorithm  returns  the  ANtVOR  graph  reinesentaticm  of  assembly  sequoices.  The  correcmess  of  the 
algorithm  is  based  cm  the  assumption  that  it  is  always  possible  to  decide  correctly  whether  two  subassemblies  can  be 
joined,  based  on  geometrical  and  physical  criteria.  liiis  paper  presents  an  approach  to  compute  this  decision.  An 
experimental  implementation  for  the  class  of  products  m^le  up  of  polyhedral  and  cylindrical  parts  having  planar  or 
cylindrical  contacts  among  themselves  is  described.  Bounds  for  the  amount  of  computation  involved  are  presented. 
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1.  Introduction 


The  choice  of  the  sequence  in  which  pans  or  subassemblies  are  put  together  in  the  mechanical  assembly  of  a 
product  can  drastically  affect  the  efficiency  of  the  assembly  process.  For  example,  one  sequence  may  require  less 
fixtuiing,  less  changing  of  tools,  and  include  simpler  and  mote  reliable  operations  than  othos.  The  choice  of  the 
assembly  sequence  is  usually  made  by  a  human  expert  In  the  case  of  manufacturing,  the  choice  is  typically  made 
by  an  industrial  engineer.  In  the  case  of  repair,  the  choice  is  made  by  the  maintenance  posonnel.  No  clear 
systematic  procedure  seems  to  be  followed  in  either  case.  Humans  seem  to  use  common  sense  and  past  expoience 
blended  in  a  fuzzy,  somedmes  inconsistent  and  not  well  understood  way. 

There  is  a  growing  need  to  systematize  and  to  computerize  the  generation  of  assembly  sequences  for  several 
reasons: 

•  Although  many  experienced  industrial  engineers  have  a  knack  for  devising  efficient  ways  to  assemble  a 
given  product  systematic  procedures  are  needed  to  guarantee  that  no  good  assembly  sequence  has  been 
overlooked.  For  complex  products,  the  number  of  feasible  assembly  sequences  may  be  so  large  that 
even  skillful  engineers  may  fail  to  notice  many  possibilities.  The  availability  of  a  systematic  procedure 
that  is  proven  correct  and  complete  will  guarantee  that  all  feasible  sequences  and  only  the  feasible 
sequences  will  be  generated. 

•  The  planning  and  inogramming  chores  in  manufacturing  are  time  consuming  and  enor-ptone.  For 
small  batches  of  production,  the  cost  of  planning  and  programming  can  weigh  heavily  in  the  total 
production  cost  Moreover,  the  time  spent  in  plaiming  and  programming  may  excessively  delay  the 
actual  production.  The  automation  of  these  chores  will  expedite  their  execution,  reduce  their  cost,  and 
improve  their  quality.  Systematic  {uocedures  are  needed  in  order  to  facilitate  the  automation  of 
plaiming  and  programming  of  assembly  systems. 

•  In  simultaneous  engineering  environments,  the  automation  of  sequence  planning  will  help  the  designer 
to  assess  the  assembly  process  requirements  of  diffoent  design  solutions  for  a  given  product  For  some 
products,  small  changes  in  the  design  can  have  a  large  impact  on  the  assembly  alternatives. 

•  Autonomous  systems  for  triplications  such  as  space  or  deep  sea  exploration  will  need  the  ability  to 
generate  assembly  or  disassonbly  sequences  that  fit  the  particular  situation  they  encounter.  It  is 
virtually  impossible  to  preprogram  all  possible  situations  those  systems  might  face,  particularly  if 
execution  erttna  can  occur  and  the  systems  are  expected  to  recover  autonomously. 

•  In  less  structured,  more  dynamic  manufacturing  systems  or  facilities  there  is  a  need  to  ad^t  the 
assembly  process  to  different  machines.  The  need  to  produce  difierent  products  in  the  same  shop  may 
lead  to  the  choice  of  an  assembly  sequence  for  a  product  that  may  not  be  the  most  efficient  but  uses  the 
idle  equipment  in  the  shop.  Knowledge  of  all  assembly  sequence  options  of  each  product  is  needed  in 
order  to  optimize  the  overall  use  of  machines  and  tools.  Similarly,  when  the  same  product  is  assembled 
in  different  shops,  the  knowledge  of  all  assembly  sequences  is  needed  in  the  selection  of  the  sequence 
more  suitable  to  the  equipment  available  in  each  shop. 

This  paper  presents  an  algorithm  for  the  generation  of  mechanical  assembly  sequences  and  a  proof  of  its 
correctness  and  completeness.  The  algorithm  takes  a  description  of  the  product  and  returns  the  corresponding 
AN0A3R  graph  rqnesentation  of  assembly  sequences  [17].  It  is  assumed  that  exactly  two  parts  or  subassemblies  are 
joined  at  each  time,  and  that  after  parts  have  been  put  together  they  remain  together.  It  is  also  assumed  that 
whenever  parts  are  joined  forming  a  subassembly,  all  contacts  between  the  parts  in  that  subassembly  are  established. 
These  assumptions  are  consistent  with  the  trend  towards  product  designs  that  are  suitable  for  automatic 
assembly  [2, 5]. 
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The  conectness  of  the  algorithm  is  based  on  the  assumption  that  it  is  always  possible  to  decide  correctly  whether 
two  subassemblies  can  be  joined,  based  on  geometrical  and  physical  criteria.  This  paper  presents  an  approach  to 
cnmpit<»  this  decision.  An  experimental  implementation  for  the  class  of  products  made  up  of  polyhedral  and 
cylindrical  parts  having  planar  or  cylindrical  contacts  among  themselves  is  described. 

The  apww't  of  computation  involved  in  generating  the  AND/OR  graph  rqrresentation  of  assembly  plans  depends 
on  the  niimlvtr  of  parts  that  make  up  the  product,  on  how  those  parts  are  interconnected,  and  also  on  the  resulting 
AND/OR  graph.  Bounds  for  the  amount  of  computation  involved  are  presented. 


2.  Background 

The  algorithm  presented  in  this  paper  takes  as  input  a  representadmi  of  the  product,  and  generates  the  set  of  all 
feasible  assembly  sequences  which  is  rqxesented  as  an  AND/OR  graph.  This  section  reviews  previous  work  on 
modeling  assemblies,  on  representing  assembly  sequences,  and  on  generating  assembly  sequences. 


2.1.  Modeling  assemblies 

The  research  on  high  level  l^mgnagea  far  robotic  assembly  has  explored  the  use  of  assembly  models.  That 
research  aimed  at  the  antnmatir  generation  of  the  actions  that  a  robot  should  poform  in  nder  to  assemble  a  product 
Typically  the  sequence  in  which  parts  should  be  put  together  was  given. 

One  of  the  earliest  works  cm  robot  programming  was  the  RAPT  [1]  system  in  which  bodies  were  described  in 
terms  of  their  features  such  as  planar  faces,  shafts,  and  holes.  The  spatial  relationships  between  parts  were  described 
by  triples  {type-^f-~spadal-relation  ,feature.l  ,feature.2).  Vor  example,  (fits  describes  the  spatial 

relationship  between  the  shaft  5,-  and  the  hole  Hj.  The  s^  of  q;)atial  relations  between  parts  was  input  to  an  inference 
engine,  and  the  relative  positions  of  parts  or  their  degrees  of  freedom  were  determined.  Later  extensions  to 
RAPT  [28. 29]  allowed  the  user  to  describe  assemblies  not  only  by  the  qntial  relationships  between  the  parts  but  also 
by  the  actions  required  to  bring  those  parts  together. 

Taylor  [37]  developed  a  rqxesentation  of  assemblies  based  on  attribute  grapiis.  The  nodes  in  these  graphs 
OHreqxxKl  to  either  objects,  or  features  of  objects.  Entities  that  have  volume  such  as  assemblies  and  parts  are 
objects,  whereas  entities  that  do  not  have  volume  such  as  surfaces  and  edges  are  features.  Each  link  in  the  graph 
associates  rare  node  either  to  another  node  or  to  a  link.  For  example,  a  subpart  link  may  associate  a  part,  which  is  an 
object  node,  to  an  assembly,  which  is  amther  object  node;  and  a  nominal-transformation  link  may  associate  a 
feature  node  containing  a  4x4  homogeneous  coordinate  transform  matrix  to  a  subpart  link.  The  information 
describing  the  shape  of  an  object  is  contained  either  in  the  node  corresponding  to  that  object,  if  the  shape  is  simple, 
or  in  the  nodes  corresponding  to  its  subparts,  if  the  shape  is  complex.  In  the  latter  case,  which  is  the  case  of 
assemblies,  the  composition  of  the  subpart’s  shapes  may  be  described  either  by  homogeneous  transform  feature 
nodes  associated  to  the  subpart  links,  or  by  associations  of  features  of  subparts  corresponding  to  spatial  relationships 
between  those  features.  Taylor  allows  redundancy  of  shape  description  and  both  types  of  descripdons  for  the 
compositions  of  shapes  may  coexist 

In  the  AUTOPASS  system  [39],  the  rqrresentation  of  asremblies  was  based  on  a  graph  structure  in  which  each  node 
represented  a  volumetric  entity,  either  a  part  or  a  sub-part  or  an  assembly,  and  the  edges  were  directed  and  labeled 
to  indicate  four  kinds  of  relationships:  part-of,  attachment  constraint  and  assembly-component  The  nodes  had 
attributes  which  included  the  volumetric  description  and  the  location  of  the  corresponding  object  The  part-of 
relationship  induced  a  tree  structure  on  the  assembly  model. 
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Unlike  the  woik  described  above,  which  aimed  at  high  level  languages  for  robotic  assembly,  the  work  of 
Bomjault  [6]  aimed  at  modeling  the  assembly  process.  Towards  that  goal  he  used  two  types  of  graphs  to  represent 
products.  The  graph  of  contacts  ("graphe  de  liaisons  mdcaniques”  [6])  contains  one  node  for  each  part  in  the 
assembly,  and  one  edge  ftn’  each  contact  between  two  parts.  Since  the  same  pair  of  parts  may  have  more  than  one 
contact,  the  graph  of  contacts  is  not  necessarily  simple.  From  the  graph  of  contacts,  Bouijault  defined  the  graph  of 
connections  ("gr^he  de  liaisons  fonctionelles"  [6])  which  has  one  node  fa*  each  part  in  the  assembly,  and  one  edge 
for  each  pair  of  parts  that  have  at  least  one  contact  By  definition,  the  graph  of  connections  is  always  a  simple 
graph. 

The  nsodel  of  assemUies  presented  in  section  3  is  similar  to  the  attributed  graph  used  previously  [37, 39]  but 
extended  to  incorporate  the  attachment  of  contacts.  This  extension  is  needed  to  make  possible  the  reasoning  about 
the  feasibility  of  assembly  tasks. 


22.  Representation  of  assembly  sequences 

One  assembly  sequence  can  be  represented  by  an  ordered  list  of  tasks;  therefore  it  is  possible  to  represent  the  set 
of  all  assembly  sequences  by  a  set  of  lists,  each  ccuresponding  to  a  different  assembly  sequence.  Since  many 
assembly  sequences  share  common  subsequences,  attempts  have  been  made  to  create  more  compact  representations. 

One  eariy  attempt  was  the  use  of  a  set  of  tasks  and  a  set  of  precedence  constraints  relating  two  task«  [is].  But  as 
discussed  elsewhere  [17],  there  are  products  for  which  standard  precedence  constraints  cannot  encompass  all 
sequences. 

Directed  graphs  of  assembly  states  can  explicitly  encompass  the  set  of  all  assembly  sequences.  The  nodes  in 
these  graphs  may  be  either  a  partition  of  the  set  of  parts  [18],  or  a  subset  of  connections  of  pairs  of  parts  [6, 1 1]. 
Figure  1  shows  a  directed  graph  of  assembly  states  for  a  three  part  product.  The  nodes  in  figure  1  are  labeled  by  the 
partitions  of  the  set  of  parts  containing  the  subsets  of  parts  of  each  subassembly  already  assembled  at  each  state  of 
the  assembly  process.  Lower  and  upper  bounds  for  the  size  of  these  graphs  as  a  function  of  the  number  of  parts  in 
the  product  are  presented  elsewhere  [18]. 

ANtVOR  [17]  graphs  of  subassemblies  can  also  encompass  the  set  of  all  assembly  sequences.  The  nodes  in  these 
AND/OR  graphs  conespond  to  subassemblies  and  the  hyperarcs  correspond  to  assembly  taAs  in  which  two 
subassemblies  are  joined  to  yield  a  larger  more  complex  subassembly.  The  hyperarcs  point  from  the  node 
corresponding  to  the  larger  subassembly  to  the  nodes  cmiesponding  to  the  smaller  subassemblies.  Figure  2  shows 
the  AND/OR  graph  of  subassemblies  for  a  three-part  product  The  nodes  in  figure  2  are  labeled  by  the  set  of  parts  that 
make  up  their  corresponding  subassemblies. 

Although  for  three-part  assemblies  the  AND/OR  graph  has  more  nodes  than  the  directed  graph  of  assembly  states, 
for  assemblies  with  large  number  of  parts  the  AND/OR  graph  has  substantially  fewer  nodes  than  the  directed  graph  of 
assembly  states.  Moreover,  the  AND/OR  graph  of  subassemblies  shows  explicitly  the  possibility  of  parallel  execution 
of  assembly  tasks.  Lower  and  upper  bounds  for  the  size  of  these  AND/OR  graphs  as  a  function  of  the  number  of  parts 
in  the  (voduct  are  presented  elsewhere  [18]. 

Bouijault  [6]  showed  that  a  set  of  logical  expressions  can  be  used  to  encode  the  directed  graph  of  assembly  states. 
For  a  product  that  has  L  coimections  between  pairs  of  parts,  Bourjault  represented  each  state  in  the  directed  graph  of 
assembly  states  by  a  binary  vector  £=  [  A.j ,  X2 ,  •  •  •  ,  in  which  the  1*  component  is  true  or  false  respectively  if  the 

connection  is  established  in  that  state  or  not  Let  5,-  be  the  set  of  states  from  which  the  connection  can  be 
established  without  precluding  the  completion  of  the  assembly.  Clearly,  if  5,-  has  ilT  elements,  each  element  satisfies 
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Figure  1:  The  directed  gnq)h  of  assembly  states  of  a  three-part  assembly 


Figure  2:  The  AND/OR  graph  for  a  three-part  assembly 


where  both  the  sum  and  the  product  are  logical  operations,  and  is  either  the  symbol  if  the  A:*''  element  in  5,- 
corresponds  to  a  state  in  which  the  connection  has  been  established,  or  the  symbol  X,  if  the  element  in 

corresponds  to  a  state  in  which  the  connection  has  not  been  established.  Bouijault  calls  the  left  side  of  the 

equation  above  the  establishment  condition  for  the  connection  ("conditions  de  r6alisabilit6"  [6])  and  he  represents 
it  by  C,.  It  is  often  possible  to  simplify  the  expression  of  C,  using  the  rules  of  boolean  algebra.  For  some  products, 
the  simplified  expressions  are  very  short  By  knowing  the  establishment  conditions  for  all  connections  of  a  product, 
as  well  as  the  product’s  graph  of  connections,  one  can  reccxistruct  the  directed  graph  of  assembly  states. 


2  Generation  of  assembly  plans 

Planning  has  been  an  important  research  topic  in  artificial  intelligence,  and  the  AI  aiqproach  has  dominated  much 
of  the  research  in  robot  task  planning  using  domain-independent  methods.  The  central  idea  of  domain  independent 
planning  is  to  have  one  general  purpose  inference  engine  which  can  be  used  for  any  domain  by  describing  the  initial 
state,  the  goal,  and  the  operators  in  a  logic  formalism.  But  domain-independent  planners  have  serious  limitations 
that  preclude  their  use  in  generating  assembly  sequences  based  on  a  description  of  the  product  Chapman  [10] 
reviews  the  literature  on  domain-independent  planning  and  discusses  their  main  limitations. 
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Bomjaiik  [6]  has  explored  ways  to  obtain  the  establishment  conditions  C,-  without  enumerating  all  the  states  in  the 
directed  graph  of  assonbly  states.  For  example,  he  noticed  that  an  affirmative  answer  to  the  question 

•  Is  it  true  that  the  connection  cannot  be  established  if  the  j***  connection  has  already  been  established? 

means  that  no  C  will  contain  expiessicms  including  Xj,  and  that  Cj  will  not  contain  expressions  including 
unless  the  and  the  connectimis  can  be  established  simultaneously.  Bourjault’s  method  uses  a  cleverly  chosen 
sequeiu:e  of  questions  that  can  in  many  cases  expedite  the  obtainment  of  the  establishment  conditions  for  all 
ctmnections. 

De  Fazio  and  Whitney  [11]  proposed  a  set  of  questiais  smaller  than  that  used  by  Bouijault  The  user  of  their 
method  must  first  draw  the  graph  of  cormections  corresponding  to  the  assembly,  and  then  answer  a  pair  of  questions 
for  each  cormection.  For  the  connection,  the  questions  are: 

1.  what  connections  must  be  done  prior  to  doing  the  ^  connection? 

2.  what  connecticxis  must  be  left  to  be  done  after  doing  the  ^  connection? 

The  answers  to  these  questions  should  be  expressed  in  the  fcrnn  of  precedence  relationships  between  connections  or 
between  logical  ccxnbinations  of  connections.  For  example,  the  answers  to  the  two  questions  fcv  the  connection 
C,- could  be: 

1.  (Cj  or  (C*  and  CJ)  ->  C,- 

2.  Ci  (C,  or  (C,  and  C„)) 

The  symbol  reads  must  precede,  and  Cj,  C„,  C,,  aitd  are  other  connections  between  parts  of  the 

assembly.  Once  the  precedence  relationships  have  been  generated,  a  computer  program  can  generate  the  assembly 
sequences.  Lui  [25]  describes  a  program  that  generates  the  assembly  sequences  based  on  the  precedence 
relationships  and  on  the  graph  of  connections. 

Both  of  these  approaches  [6, 11]  lend  themselves  to  interactive  systems  in  which  a  computer  profnm  generates 
the  questions,  a  human  expert  supplies  the  answers,  and  the  program  then  generates  the  precedence  relationships 
between  connections  or  between  logical  combinations  of  connections.  For  simple  cases,  these  approaches  have  the 
advantage  that  they  exploit  the  engineer’s  intuitive  understanding  of  parts  relations  and  feasibility  of  operations.  For 
complex  cases,  it  may  be  very  difficult  for  a  human  expert  to  answer  the  questions  and  to  guarantee  the  cotiecmess 
of  the  answers.  And  even  assuming  that  the  qiMStions  are  answered  correctly,  proofs  ot  cotiecmess  and 
completeness  of  the  algorithms  are  needed  to  guarantee  that  the  resulting  jnecedence  relations  are  satisfied  by  all  the 
feasible  assembly  sequences  and  only  by  the  feasible  assembly  sequences.  Neither  Bourjault  nor  De  Fazio  and 
Whimey  have  formally  proven  the  cotiecmess  and  completeness  of  their  algtnthms. 

Furthermore,  it  seems  very  difficult  to  develop  computer  programs  that  will  answer  the  questions  in  eitho’  method 
directly  from  a  description  of  the  assembly;  any  system  based  on  these  tqiproaches  will  need  the  human  expert  to 
supply  the  answers.  In  the  cases  in  which  precedence  relationships,  together  with  the  assembly’s  graph  of 
connections  provide  a  useful  representation  of  assembly  sequences,  an  alternative  to  have  the  questions  answered  by 
a  human  expat  is  to  have  them  answered  by  a  program  that  takes  as  input  the  set  of  assembly  sequences  generated 
by  the  algorithm  presented  in  this  paper. 

3.  A  Relational  Model  for  Assemblies 

A  mechanical  assembly  is  a  composition  of  parts  intoconnected  forming  a  stable  unit  Each  part  is  a  solid  object 
Parts  are  interconnect^'!  whenever  they  have  one  or  more  surfaces  in  contact  Surface  contacts  between  parts  reduce 
the  degrees  of  freedom  for  relative  motion.  A  cylindrical  contact  for  example,  prevents  any  relative  motion  that  is 
not  a  translation  along  the  axis  or  a  rotation  around  the  axis.  Attachments  may  act  on  surface  contacts  and  eliminate 
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all  degrees  of  freedom  for  relative  motion.  For  example,  if  a  cylindrical  contact  has  a  pressure-fit  attachment,  then 
no  relative  motion  between  the  parts  is  possible. 

The  representations  of  products  developed  for  high  level  robot  programming  languages  emphasized  the  geometric 
aspects  such  as  the  shape  of  the  parts  and  the  contacts  between  parts.  That  emphasis  is  consistent  with  the  goal  of 
generating  a  sequence  of  robot  actions  that  will  join  two  subassemblies,  given  the  sequence  in  which  parts  or 
subassemblies  should  be  put  together.  However  for  the  generation  of  the  assembly  sequences,  a  purely  geometric 
description  of  the  produa  is  not  sufficient  There  are  sequences  that  would  be  feasible  frx»n  a  geometric  point  of 
view,  but  are  unfeasible  in  practice  due  to  forces  resulting  from  fasteners.  Therefore,  a  model  of  assemblies  to  be 
used  in  generating  assembly  sequences  must  represent  explicitly  the  fastenings  that  bind  one  part  to  another. 

The  representation  of  assemblies  used  by  the  algorithms  described  in  sections  4  and  S  is  a  relational  model  that 
includes  three  types  of  entities:  parts,  contacts,  and  attachments.  It  also  includes  a  set  of  relationships  between 
entities.  Both  entities  and  relationships  can  have  attributes.  Fbrmally,  the  relational  model  of  an  assembly  is  a 
S-tuple  {P  ,C  ,  A  ,R  ,  (t-functions)  in  which 

•  F  is  a  set  of  symbols,  each  of  which  corresponds  to  one  part  in  the  assembly.  No  two  elements  of  P 
correspond  to  the  same  part 

•  C  is  a  set  of  symbols,  each  of  which  corresponds  to  a  contact  between  surfaces  of  two  parts  of  the 
assembly.  No  two  elements  of  C  correspond  to  the  same  ccmtact  The  two  surfaces  must  be  compatible. 

'  An  example  of  a  pair  of  compatible  surfaces  are  a  cylindrical  shaft  and  a  cylindrical  hole.  The  same 
pair  of  parts  may  have  more  than  one  contact  And  the  same  surface  of  one  part  may  be  in  contact  with 
surfaces  of  two  or  more  other  parts. 

•  A  is  a  set  of  symbols,  each  of  which  corresponds  to  an  attachment  that  acts  on  a  set  of  contacts.  No  two 
elements  of  A  correspond  to  the  same  attachment  An  attachment  always  has  an  agent  which  can  be 
either  the  attached  contact  or  another  contact  or  a  part  The  access  to  an  attachment  may  be  blocked 
by  one  or  more  parts. 

•  R  is  a  set  of  symbob,  each  of  which  corre^nds  to  a  relationship  between  pairs  of  elements  of 
PkjCuA.  No  two  elements  off?  correspond  to  the  same  relationship. 

•  a-functions  b  a  set  of  attribute  functions^  whose  domains  are  subsets  of  FuCuAu/?.  These 
functions  associate  entities  or  relationships  to  their  characteristics  such  as  the  type  of  attachment,  the 
entities  related  by  a  relationship,  and  the  stuqw  of  a  part 

This  definition  of  a  relational  model  representation  of  assemblies  is  sufficiently  general  to  encompass  a  large  class 
of  assemblies.  The  set  of  functions  can  be  enlarged  to  include  all  the  information  that  might  be  necessary  to 
generate  assembly  sequences.  In  practice,  it  may  be  convenient  to  restrict  the  class  of  assemblies  represented.  Our 
current  expoim^tal  implementation  has  the  following  restrictions: 

•  The  contacb  between  parts  involve  one  of  the  following  pairs  of  compatible  surfaces: 

•  planar  surface  and  another  planar  surface, 

*  cylindrical  shaft  and  cylindrical  hole, 

*  polyhedral  shaft  and  polyhedral  hole, 

•  threaded  cylindrical  shaft  and  threaded  cylindrical  hole. 


fimction  is  defined  is  •  subset  of  the  caitesian  product  of  two  seu  (the  domiin  tnd  the  nnge)  that  has  no  two  pain  whose  Tint  elements  ate 
the  same,  and  such  that  every  element  in  the  domain  appean  in  one  pair. 
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The  types  of  attachments  are: 

•  glue  attachment, 

•  pressure  fit  attachment, 

•  clip  attachment, 

•  screw  attachment 

The  attribute  functions  are  the  following: 

•  The  function  that  associates  a  part  to  a  description  of  its  shape: 

shapeiP  ->  S 

where  5  is  the  set  of  all  shape  descriptions. 

•  The  function  that  associates  a  part  to  a  description  of  its  location: 

location’.P  T 

where  ristheset(^all4x4  homogeneous  transformation  matrices.  The  matrix  T,-  associated  to 
part  Pf,  corresponds  to  the  position  and  orientation  of  a  reference  frame  attached  to  part  p,-  with 
respect  to  a  global  frame  of  reference  for  the  whole  assembly.  The  choice  of  this  global  frame  of 
reference  is  arbiuary,  but  the  same  global  reference  must  be  used  for  all  parts. 

•  The  fimction  that  associates  a  contact  to  its  type: 

type-of-cotuact:C  contact-types 

where  contact-types^  ( planar  ,  cylindrical ,  slot ,  threaded-cylindrical }. 

•  The  function  that  associates  a  planar  contact  to  the  coordinates,  with  respect  to  the  assembly’s 
global  frame  of  reference,  of  a  vector  normal  to  the  contact  plane 

normal'.  { c|  [c  e  C]  a  [type~of-contact{c)-p]aieit] )  -> 

•  The  function  that  associates  a  planar  contact  to  the  part-relationship  that  relates  the  contact  to  the 
part  that  is  back  of  the  contact: 

back:  { c|  [c  6  C]  a  [rype-o/-co/iracr(c)= planar] )  R 
This  function  must  be  consistent  with  the  function  normal. 

•  The  function  that  associates  a  planar  contact  to  the  part-relationship  that  relates  the  contact  to  the 
part  that  is  forward  of  the  contact 

forward:  { c|  [c  e  C]  a  [(yp«-r^-conract(c)= planar] }  -> R 
This  fimction  must  be  consistent  with  the  function  normal. 

•  The  function  that  associates  a  cylindrical,  slot,  or  threaded-cylindrical  contact  to  the  coordinates, 
with  respect  to  the  assembly’s  global  frame  of  reference,  of  the  line  of  the  axis  of  both  the  hole 
and  the  shaft 

axis:  {c  |  [c  e  C]  a 

[  type-qf-contact  ( c )  e  { cylindrical ,  slot ,  threaded-cylindrical )  ] )  -»  x 

•  The  function  that  associates  an  attachment  to  its  type: 

type-of-attachment:A  -*  attachment-types 
where  attachment-types^  ( clip  ,  pressure  ,  screw ,  glue }. 

•  The  function  that  associates  a  relationship  to  its  type: 

type-of-relationship:R  -*  relationship-types 

where  relationship-types  =  {  part-contact  target-attachment  agent-attachment  blocking-part- 
attachment  } 


•  The  fiinf-rirm  that  gg«irintftg  a  part  or  a  contact  to  its  part-contact  relationships: 

part-contact-relationships  :PkjC  ->  n(/?) 
where  n  (/{ )  is  the  set  of  all  subsets  of  R. 

•  The  function  that  associates  a  part-contact  relationship  to  its  part: 

part:[r\[r^  /?]  a  [ope-o^-relflrio#iiAi>(r)=part-contact] }  P 

•  The  function  that  associates  a  part-contact  relationship  to  its  contact: 

contact:  { rl  [r  e  i?]  a  [ope-no/-r«tot*o»jWp(r)= part-contact] )  -»  C 

•  The  finvfifw  that  agwiattiia  an  attachment  or  a  contact  to  its  target-attachment  relationships: 

target-attachment-relationships:  CkjA  ->  n(J?) 

•The  function  that  associates  an  attachment,  a  contact  or  a  part  to  its  agent-attachment 
relationships: 

agent-attachment-relationships'.PKjCKjA  ->!!(/?) 

•  The  fimction  that  associates  a  taiget-attachmeht  relationship  to  its  contact 

target:  { r|  [r  e  /?]  a  [type-of-relationship(r)=tatgtt-atsacittaexA]  ]  -*  C 

•  The  function  that  associates  an  agent-attachment  relatimiship  to  its  agent 

agent:  { r]  [r  6  1?]  a  [type-of-relationshipirysAgent-attacbment]  ]  -*  PuC 

•  The  fimetiQn  that  a  bhxking-^yart-attaclunent  relationsh^  to  its  blocking-part 

Nocking-part:  [r\[r  e  i?]  a 

A.[type-of-relationship(r)»blocldag-pm~maciaoent]]  P 
•The  function  that  assodates  a  target-attachment,  a  blocking-part-attachment,  or  an  agent- 
attachment  lelatitMish^  to  its  attachment 

■  attachment:  {r|tr  €  i?]  a  [type-<^-relationship{r)  6  J5] )  ->  A 
with^=(  target-attachment, blocking-part-attachment, agent-attachment  ). 

The  relational  model  of  an  assembly  must  be  consistent  For  example,  if  part(rj  )spj  and  contact  {,r^)sc^  that 
r^e  part-contact-relationships {pi)  and  r^^^  part-contact-relationships (p<^)  must  hold.  Furthermore,  the 
relational  model  of  an  assembly  must  satisfy  some  syntactic  constraints,  the  most  important  of  which  are: 

•  every  contact  most  have  exactly  two  part-contact  relationships; 

•  every  part  most  have  at  least  one  part-contact  relationshq),  excejpt  in  the  case  the  assembly  has  only  one 
part; 

•  every  attachment  must  have  at  least  one  target-attachment  relationship,  and  at  least  one  agent- 
attachment  relationship. 

The  relational  model  of  an  assembly  can  be  rqiresented  by  a  gr^h  plus  the  associated  attribute  functions.  Figure 
3  shows  a  simple  product,  and  figure  4  shows  its  anresponding  relational  model  graph. 

The  nodes  in  figure  4  correspond  to  the  entities.  Nodes  corresponding  to  part  entities  are  rectangles,  nodes 
corresponding  to  contact  entities  are  circles,  and  nodes  corresponding  to  attachment  entities  are  triangles.  All  nodes 
contain  labels  indicating  their  cmresponding  entities.  The  attribute  functions  associated  with  the  contact  entities  are 
shown  in  Table  1. 

The  labeled  lii.  s  connecting  two  nodes  in  figure  4  correspond  to  the  relationships.  Excqtt  for  R5,  R6,  R13,  and 
R14,  all  relationships  are  part-contact  Relationships  R5  and  R13  are  target-attachment;  they  indicate  that  the 
contacts  C2  and  C5,  respectively,  are  attached.  Relationships  R6  and  R14  are  agent-attachment;  they  indicate  that  the 
agents  of  the  attachments  are  the  target  contacts  themselves.  Next  section  (see  figures  9  and  10)  shows  an  example 
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CAP 


STICK 


RECEPTACLE 


HANDLE 


Figures:  A  simple  product  in  exploded  view 


Figure  4:  The  relational  model  graph  for  the  product  show  in  figure  3 
of  an  attachment  whose  agent  is  not  the  contact  itself. 

Given  the  relational  model  of  a  product  {P  ,C  ,  A  ,R ,  a-fimcAons),  a  number  of  otho’  useful  rqnesentations 
may  be  generated.  For  example,  the  graph  of  connections  of  the  assembly,  as  defined  by  Bouijault  [6]  (see  section 
2),  is  the  simple  graph  {V,E)'m  which 

V=P 

E  =  {  (P, ■./»/)  I  \Pi  6  P]  A  \pj  e  P]  A 

aBcBtj  3r2  [  [c  e  C]  A  [{rj,r2}=/>arr-cortract-retotK?/«Aips(c)]  a 

A  [p,=part(r,)]  a  [pj=part(r2)'\  ]  ) 

Figure  S  shows  the  graph  of  connections  for  the  simple  product  shown  in  figure  3. 
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Table  1:  Attribute  Functions  for  the  Contact  Entities  in  Figure  4 


Cl 

C2 

C3 

V 

C5 

type-of-contact 

planar 

threaded- 

cylindrical 

cylindrical 

planar 

threaded- 

cylindrical 

normal 

(010) 

nil 

nil 

(010) 

nil 

back 

cap 

nil 

nil 

stick 

nil 

forward 

stick 

nil 

nil 

handle 

nil 

axis 

nil 

((000X010)) 

((000)  (010)) 

nil 

((000)  (010)) 

part-contact 

relationships 

(R1  R2) 

(R3R4) 

(R7R8) 

(RO  RIO) 

(R11  R12) 

target-attachments 

relationships 

nil 

(R5) 

nil 

nil 

(R13) 

agent-attachment 

relationships 

nil 

(R6) 

nil 

nil 

(R14) 

C  s  cap  S  s  stick  R  » lecqTtacle  H  s  handle 
Figure  5:  The  gisq)fa  of  connections  for  the  product  shown  in  Hgure  3 


3.1.  Subassemblies 

A  subassembly  is  a  nonempty  subset  of  parts  that  either  has  only  one  element  (i.e.  only  (xie  part),  or  is  such  that 
every  part  has  at  least  one  surface  contact  with  another  part  in  the  subset  Although  there  are  cases  in  which  it  is 
possible  to  join  the  same  pair  of  parts  in  more  than  one  way,  a  unique  assembly  geometry  will  be  assumed  for  each 
pair  of  parts.  This  geometry  corresponds  to  their  relative  location  in  the  whole  assembly.  A  subassembly  is  said  to 
be  stable  if  its  parts  maintain  their  relative  position  and  do  not  break  contact  spontaneously.  All  one-part 
subassemblies  are  stable. 

Given  the  relational  model  of  a  product  {P  ,C ,  A  ,R  ,  a-functions),  the  relational  model  of  a  subassembly  of 
that  product  is  a  relational  model  {Pg  .Cg  ,Ag  ,Rg  ,  a-functions^)  in  which  P^  c  F,  c  C,  ^A,R^^  R, 
and  every  function  in  a-functions^,  is  a  subset  of  the  corresponding  function  in  a-functions.  In  addition  to  the 
syntactic  constraints  mentioned  above  that  every  relational  model  of  an  assembly  must  satisfy,  the  relational  model 
{Pg  ,Cs  ,A^  ,Rg  ,  a-functions f)  of  a  subassembly  oi{P  ,C  ,A  ,R  ,  a-functions)  must  also  satisfy  the  constraint: 
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Vc  Vtj  Vtj  [  [c  e  C]  A  [  { Tj  .Tj }  =‘part-coniact-relationships(,c)]  a 

A  [pan{r{)e  a  [partir^)  €  F^]  ]  -»  [c  e  C^] 

This  constraint  corresponds  to  the  assumption  that  wh^ver  parts  are  joined  forming  a  subassembly  all  contacts 
between  the  parts  in  that  subassembly  are  established.  It  requires  that  those  contacts  in  the  model  of  the  assembly 
whose  two  part-contact  relationships  involve  parts  in  the  subassembly  must  also  be  in  the  model  of  the  subassembly. 
For  example,  to  the  product  shown  in  figure  3,  there  is  no  subassembly  relational  model  in  which  {  CAP. 
RECEPTACLE,  STICK  ],  and  C^-  {  C2,  C3  }.  If  both  the  cap  and  the  stick  are  in  F^,  then  contact  Cl  most  also  be  in 
Cg.  This  constraint  allow  the  characterization  of  any  subassembly  <F^  .Cg  ,Ag  ,Rg  ,  a-functionsg)  of  a  product 
{P  ,C  ,A  ,R  ,  a-fioictions)  by  its  set  of  parts  Pg  only.  This  feanue  will  be  used  in  the  algorithm  to  the  generation 
of  mechanical  assembly  sequences  described  in  the  subsequent  sections.  In  that  algorithm,  the  intermediate 
subassemblies  will  be  characterized  by  their  sets  of  parts.  Given  a  subset  of  parts  Pg,  there  is  a  corresponding 
subgnqph  <  Vg  ,Eg)  oi  the  assembly’s  grtq)h  of  connections  (V,E).  In  this  subgraph,  the  set  of  nodes  Vg  includes  all 
the  elements  of  V  that  ccsrespond  to  the  parts  in  V^.  And  the  set  of  edges  Eg  includes  all  the  elements  of  E  that  have 
both  end  points  in  Vg.  A  subset  of  parts  Pg  characterize  a  subassembly  if  and  only  if  the  conesponding  subgr^h 
<  Vg,Eg)  is  connected  (Le.  has  only  one  component).  A  predicate  that  is  satisfied  only  by  the  subsets  of  parts  that 
correspond  to  subassemblies  can  be  defined  as  follows: 

Definition  1:  The  subassembly  predicate  associated  to  subassemblies  of  assembly 
^=<F  ,C  ,A,R  ,  a-funcUons)  is  the  predicate 
sa^:  n(F)  -»  {true, false) 
with  sa^(6)=true  if  the  subgr2q)h  <  V^,£^)  in  which 

E  =  {  (.Pi.Pj)  I  iFi  6  Fj]  A  \pj  €  Fj]  A 

A  3c  3rj  3r2  [  [c  6  C]  a  [  { fj ,r2}  =part-contact-relcuionships(c)]  a 

A  [Pi=part(.r^)]  a  [pj-partir^)]  ]  } 

is  connected. 

4.  Decompositions  of  a  Relational  Model  of  an  Assembly 

The  problem  of  generating  the  assembly  sequeiKes  to  a  product  can  be  transformed  into  the  problem  of 
generating  the  disassembly  sequences  for  the  same  product  Since  assmbly  tasks  are  not  necessarily  reversible,  the 
equivalence  of  the  two  problems  will  hold  only  if  each  task  used  in  disassembly  is  the  reverse  of  a  feasible  assembly 
task,  regardless  of  whether  this  reverse  task  itself  is  feasible  or  not  The  expression  disassembly  task,  therefore, 
refers  to  the  reverse  of  a  feasible  assembly  task, 

As  mentioned  in  the  introduction,  it  was  assumed  that  exactly  two  parts  or  subassemblies  are  Joined  at  each  time. 
It  was  also  assumed  that  whenever  parts  are  joined  forming  a  subassembly,  all  contacts  between  the  parts  in  that 
subassembly  are  established.  In  the  disassembly  problem,  each  task  splits  one  subassembly  into  two  smaller 
subassemblies,  maintaining  all  contacts  between  the  parts  ui  either  of  the  smaller  subassemblies. 

A  decomposition  tqjproach  can  be  used  to  solve  the  disassembly  problem.  In  this  approach  the  problem  of 
disassembling  one  assembly  is  decomposed  into  two  distinct  subproblems,  each  being  to  disassemble  one 
subassembly.  Every  decomposition  must  correspond  to  a  disassembly  task.  If  solutions  for  both  subproblems  that 
result  fiom  the  decompositions  are  found,  a  solution  for  the  original  problem  can  then  be  obtained  by  combining  the 
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solutions  to  the  two  subinobleins  and  the  task  conesptMiding  to  the  decomposition.  For  subassemblies  that  contain 
one  part  only,  a  trivial  solution  containing  no  assembly  task  always  exists.  This  decomposition  approach  lends  itself 
to  an  AND/OR  graph  rqjresentatkm  of  assembly  sequences  [17].  The  correspondence  between  the  AND/OR  graph  and 
the  directed  graph  rqjresentations  of  assembly  sequences  is  discussed  elsewhere  [18]. 

From  now  on,  references  to  products,  to  assemblies,  or  to  subassemblies  are  references  to  their  relational  models, 
which  are  always  assumed  to  be  consistent  and  to  satisfy  the  syntactic  constraints  of  a  relational  model  of  an 
assembly.  A  teal  product  will  be  referred  to  as  a  physical  product,  a  real  assembly  as  a  physical  assembly,  and  a  real 
subassembly  as  a  physical  subassembly. 

A  decomposition  of  an  assembly  {P  ,C  ,A  ,R  ,a-fimctions)  is  a  pair  of  its  subassemblies 
{P^i  ,Csi  ,Agi  ,Rgi  ,o-fiuictionsgi)  and  <Pj2  <^52  • '^52*^52  such  that  Psi'^^si^^ 

Pgi  Th®  set  C5j_52=C-(Cjj  uC^)  is  referred  to  as  the  contacts  of  the  decomposition;  they  are  the 

contacts  that  belong  to  C  and  do  not  belong  to  either  ^51  or  Cj2-  The  contacts  of  a  decomposition  of  an  assembly 
define  a  cut-set  in  that  assembly’s  graph  of  connections.  Conversely,  a  cut-set  in  the  gr^h  of  connections  of  an 
assembly  define  a  decomposition  of  that  assembly. 

A  decomposition  of  an  assembly  is  said  to  be  feasible  if  it  satisfies  three  predicates:  GEOMETRIC-FEASIBILITY, 
MECHANICAL-FEASIBILITY,  and  STABILITY.  These  predicates  reflect  the  feasibility  of  joining  the  physical 
subassemblies  to  produce  the  physical  assembly. 

The  GEOMETRIC-FEASIBILITY  predicate  is  true  if  there  is  a  coUision-Ctee  path  to  bring  the  two  subassemblies 
into  contact  firom  a  situation  in  which  they  are  sufficimdy  far  apart  Foe  the  assembly  shown  in  figure  3,  for 
example,  there  is  no  colliskm-fiee  path  that  will  bring  the  stick  into  contact  with  the  subassembly  made  up  of  the 
cap,  the  receptacle,  and  the  handle.  Joining  the  stick  to  die  subassembly  made  up  of  the  three  other  parts  is  said  to 
be  geometrically  unfeasMe.  Joining  the  stick  to  the  subassembly  made  up  of  the  cap  and  the  receptacle,  howevo-,  is 
geometrically  feasible  since  there  is  a  coUision-firee  path  to  Ining  the  two  subassemblies  into  contact 

The  MECHANICAL-FEASIBILITY  predicate  is  true  if  it  is  feasible  to  establish  the  fftTffrhmfints  that  act  on  the 
contacts  of  the  decomposidon.  Hgure  6  shows  a  three-part  assembly  in  which  the  part  in  the  center  (part  B)  is 
attached  to  the  part  in  the  right  (part  C)  through  two  built-in  bolts.  Althou^  it  is  geometrically  feasible  to  join  the 
part  in  the  right  (part  C)  to  the  subassembly  made  iqi  of  the  two  other  parts,  it  is  impossible  to  establish  the 
attachmoits  because  the  access  to  the  bolts  is  blocked  hy  the  part  in  die  left  (part  A).  Joining  the  part  in  the  right 
(part  C)  to  the  subassembly  made  iqi  of  the  two  other  parts  is  said  to  be  mechanically  utfeasible. 

The  STABILITY  predicate  is  true  if  the  parts  in  either  physical  subassembly  maintain  their  relative  position  and  do 
not  break  contact  spontaneously.  For  the  assembly  shown  in  figure  7,  the  subassembly  made  up  of  the  parts  B  and  C 
is  not  stable  since  the  two  parts  will  break  contact  spontaneously  due  to  gravity,  regardless  of  their  cxientadon  in 
space^. 

In  practice,  the  feasibility  of  joining  two  subassemblies  depends  on  the  availability  of  adequate  resources  such  as 
machines,  tools,  and  fixtures.  For  the  general  analysis  presented  here,  it  is  assumed  that  all  such  resources  are 
available. 

As  discussed  in  secdon  3,  the  subassemblies  of  a  given  assembly  'V^{P  ,C  ,A,R ,  a-functions)  can  be 


^Ahhough  ooittactt  like  (hat  between  parts  B  and  c  in  ftgaie  7  are  not  handled  by  oor  experimental  implementation,  they  ilhiitrate  veiy  clearty 
the  stability  or  unstabOity  of  sabusembliet. 
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A  B  C 

Figure  6:  An  assmbly  that  illustrates  the  mechanical  feasibility  predicate 


Figure  7:  An  assembly  that  illustrates  the  stability  predicate 


characterized  by  their  sets  of  parts.  Therefore,  the  three  predicates  described  above  can  be  defined  as  follows: 

Definition  2:  The  geometric*feasibility  predicate  associated  to  subassemblies  of  assembly 
,C  ,  A  ,R  ,  (i-/iincdo/w>,in  which/*={pj,p2.  •  ,p^),  is  the  predicate 
n(/’)  X  n(F)  -*  (true, false} 

with  gf^(0j,62)»true  if  and  otdy  if  9jn62=0  and  diere  is  a  collision-fiee  path  to  bring  the  two 
physical  subassemblies  of  Y  characterized  by  and  82  contact  from  a  situation  in  which  they  are 
sufficiently  far  apart 

Definition  3:  The  mechanical«fearibiUty  predicate  associated  to  subassemblies  of  assembly 
'¥-{P  ,C  ,A  ,R  .a-:/i(Rcdbns),  in  which Pi  ,P2>  .Pjyl.istheiHedicate 
n(/*)  X  nCF)  -♦  (true, false) 

with  ,62)- true  if  and  only  if  9j  062=0,  ^  feasible  to  establish  the  attachments  that  act  on 

the  set  of  contacts  between  parts  in  9j  and  parts  in  62. 

Definition  4:  The  stability  predicate  associated  to  subassemblies  of  assembly 
*F=</’ ,  C  ,  A  ,  ,  a-fmctions),  in  which  /*=  (pj  ,P2»  •  *  ’  ^Pn  ) .  the  predicate 

(true, false) 

with  st^(9)sxrat  if  and  only  if  the  parts  in  6  maintain  their  relative  position  and  do  not  break  contact 
spontaneously. 

The  GEOMETRIC-FEASIBILITY  predicate  can  be  computed  using  path  planning  algorithms  [13, 20, 38)  to 
generate  a  collision-free  path  to  taing  the  two  subassemblies  into  contact,  or,  equivalently,  a  collision-fiee  path  to 
separate  the  two  subassemblies.  These  algorithms  typically  involve  large  amounts  of  computation  and  more 
efficient  approaches  to  general  path  feasibility  tests  are  needed.  For  many  industrial  assemblies,  the  computation  of 
geometric  feasibility  can  be  significantly  reduced  by  poforming  a  simple  local  analysis  which  can  indicate  that  a 
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collision>free  path  does  not  exist  For  a  given  decomposition,  this  local  analysis  looks  at  the  assembled  assembly 
and  checks  whether  there  exists  an  incremental  translation  of  one  of  the  two  subassemblies  that  is  not  blocked  by 
any  of  the  contacts  between  one  of  its  parts  and  one  of  the  parts  in  the  other  subassembly. 

Fix'  many  types  of  contacts  there  are  very  few  feasible  motions  between  the  parts.  For  example,  the  only  direction 
almig  which  a  pin  in  a  h(de  can  translate  is  the  direction  of  the  axis.  Whenever  the  part  or  subassembly  has  such  a 
constraining  contact,  the  local  analysis  can  be  performed  by  checking  the  compatibility  of  the  most  restrictive 
contact  widi  all  other  crmtacts.  In  the  case  of  the  pin  in  tire  hole,  the  local  analysis  consists  of  checking  whether  a 
translation  of  die  pin  along  its  axis  is  not  blocked  by  any  of  the  other  contacts  between  the  pin  and  the  other  part  or 
subassembly. 


The  local  analysis  is  more  difHcult  when  the  part  or  subassembly  to  be  disassembled  is  constrained  by  planar 
contacts  only.  Each  planar  contact  leaves  an  infinite  number  of  unconstrained  directions  along  which  translation  is 
possible.  All  these  directions  have  positive  [Kojection  over  the  normal  to  the  surface  of  the  blocking  part,  pointing 
towards  the  outside  of  the  blocking  part 


In  order  to  decide  whether  a  set  of  planar  contacts  does  not  completely  constrain  one  part  or  subassembly,  one 
must  find  whether  diere  is  a  nonzero  solution  to  the  system  of  linear  inequalities 
3 


X, ::  0 
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where  n  [iij  j  nj2  r,-3  ]  is  the  normal  to  the  surface  of  the  i***  contact  This  system  of  linear  inequalities  defines  a 
polyhedral  convex  cone.  It  has  been  shown  [16]  that  such  a  polyhedral  convex  cone  can  be  built  up  from  its 
(unique)  d-dimensional  face  and  its  (d  l)-dimensiofud  faces  (if  any),  where  dB3-nuik(Af ),  and  Af  is  the  matrix 
of  the  coefficients  /ijy .  If  d  is  greater  than  zero,  then  the  pol^iedral  convex  cone  has  a  face  of  dimension  greater 
than  zero  and  therefore  the  system  of  inequalities  has  a  nonzero  solution.  If  d  is  equal  to  zero,  then  the  system  of 
inequalities  has  a  nonzero  solution  only  if  the  polyhedral  crxivex  cone  has  at  least  one  one-dimensional  face.  The 
existence  of  a  one-dimensional  face  can  be  determined  by  checking  die  N’{N- 1 )  pairwise  intersections  of  the 
planes  corresponding  to  the  inequalities.  Each  intersection  of  two  distinct  planes  is  a  line.  If  one  of  the  two  unity 
vectcxs,  ^  and  -£  along  the  intersection  line  of  two  planes  has  positive  projection  over  all  the  normal  vectors 
Rj ,  /i2.  ■ '  ■  ,  rifi,  then  the  half-line  defined  by  that  vector  (£  or  -£)  is  a  one-dimensional  &ce  of  the  polyhedral 
convex  cone. 


If  thoe  is  a  nonzero  solution  to  the  system  of  inequalities,  then  the  part  (x  subassembly  is  not  completely 
constrained.  Otherwise  the  subassembly  is  completely  constrained,  and  there  is  no  need  to  look  for  a  collision  free 
path. 

Our  current  implementation  is  not  limited  to  only  checking  the  existence  of  a  nonzero  solution  to  the  system  of 
linear  inequalities,  but  includes  the  computation  of  the  polyhedral  convex  cone  of  all  solutions.  Appendix  I 
addresses  the  reasoning  about  the  feasibility  of  local  translations  for  robotic  assembly  of  a  part  constrained  by  platuir 
contacts. 


The  usefulness  of  the  local  analysis  is  further  enhanced  by  the  use  of  virtual  contacts  to  describe  blocking 
relationships  equivalent  to  contacts.  In  the  product  shown  in  figure  3,  for  example,  if  the  stick  did  not  touch  the 
handle,  the  local  analysis  as  described  above  would  indicate  that  the  stick  can  translate  (incrementally)  along  its 
axis.  In  a  case  like  this,  a  virtual  planar  contact,  analogous  to  C4  on  figure  4,  would  be  added  to  the  relational  model 
indicating  the  blocking  of  the  stick  by  the  handle. 

The  MECHANICAL-FEASIBILITY  predicate  can  be  computed  by  inspection  of  the  relational  model  of  the 
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Figure  8:  The  relational  model  the  assembly  shown  in  figure  6 

assembly.  In  our  current  implementation,  a  procedure  MECHANICAL-FEASIBILITY  checks  whether  the 
attachments  acting  on  the  contacts  of  the  decomposition  are  not  blocked  in  the  resulting  assembly,  and  are  not 
present  in  either  one  of  the  subassemUies.  Two  examples  will  illustrate  this  computation. 

Hgure  8  shows  the  relational  model  of  the  assembly  shown  in  figure  6.  Relationships  R10  and  R13  are  agent- 
attachment,  relationships  R11  and  R14  are  target-attachment,  and  relationships  R9  and  R12  are  bloddng-pait- 
attachment;  all  other  relationships  are  part-contact  The  relationships  R9  and  R12  indicate  that  part  A  blocks  the 
access  to  attachment  A2  and  to  attachment  A4.  One  of  the  disassembly  tasks  whose  mechanical  feasibility  must  be 
computed  is  the  separation  of  part  C  from  the  subassembly  made  up  of  parts  A  and  B.  The  mechanical  unfeasibility 
of  this  task,  can  be  detected  by  inspecticHi  of  the  relational  model  which  indicates  that  the  attachments  acting  on  the 
contacts  of  the  decomposition  are  blocked  by  part  A.  After  part  A  is  removed,  those  attachments  will  no  longer  be 
blocked  and  part  C  can  be  sqiarated  from  part  B. 

Figure  9  shows  an  assembly  that  has  three  parts:  a  box,  a  cover,  and  a  clip  that  attaches  the  cover  to  the  box. 
Figure  10  shows  this  assembly’s  relational  model  Relationships  R7,  R8,  and  R9  in  figure  10  are  target-attachment; 
they  indicate  that  the  three  contacts  Cl,  C2,  and  C3  are  attached  by  attachment  A1.  Relationship  R10  is  agent- 
attachment;  it  shows  that  the  agent  of  attachment  Al  is  the  clip.  One  of  disassembly  tasks  whose  mechanical 
feasibility  must  be  computed  is  the  separation  of  the  cover  fix}m  the  subassembly  made  up  of  the  box  and  the  clip. 
The  mechanical  unfeasibility  of  this  task  can  be  detected  by  inspection  of  the  relational  model  which  shows  that  the 
contacts  cannot  be  detached  while  the  agent  of  the  attachment  is  present  The  separation  of  the  clip  from  the 
subassembly  made  up  of  the  box  and  the  cover,  however,  is  feasible  because  the  agent  of  the  attachments  is  being 
separated. 

The  computation  of  the  STABILITY  predicate  will  depend  on  additional  assumptions  about  the  assembly  process. 
For  example,  it  may  be  assumed  that  all  subassemblies  can  be  made  stable  through  the  use  of  jigs  and  fixtures.  In 
our  current  implementation  we  made  this  assumption  and  we  do  not  compute  the  STABILITY  predicate.  In  previous 
work  aimed  at  selecting  an  assembly  sequence  [17],  we  assessed  the  stability  of  a  subassembly  by  the  degrees  of 
freedom  for  relative  motion  between  parts.  Similarly,  one  can  establish  a  threshold  on  the  degrees  of  freedom  for 
relative  motion  above  which  a  subassembly  would  be  considered  unstable. 
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Figure  9:  Assembly  example 


Figure  10:  Relational  model  fo’  the  assembly  example  show  in  figure  9 

An  alternative  approach  for  the  computation  of  the  STABILITY  predicate  is  to  check  whether  there  is  an 
orientation  of  the  subassembly  such  that  there  is  no  relative  motion  between  parts  due  to  gravity.  As  a  first 
rqrproximation,  friction  can  be  ignrxed  since  it  typically  helps  the  stability.  Boneschanscher  et  al.  [4]  have  taken  this 
approach  with  the  additional  assumption  that  the  subassembly  sits  on  a  table.  They  used  a  convex  hull  algorithm  to 
find  candidate  orientations  in  which  the  subassembly  can  sit  on  a  table,  and  for  these  orientations  they  checked  the 
static  stability.  Their  analysis  takes  friction  into  account 

For  the  discussion  in  the  next  section,  which  presents  the  algorithm  for  generating  the  assembly  sequences,  it  is 
assumed  that  there  exist  correct  algorithms  for  computing  the  three  predicates  discussed  above,  and  that  they  are 
combined  into  the  procedure  FEASIBIUTY-TEST  shown  in  figure  1 1. 
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procedure  FEASIBIUTY-TEST{decomposition,  assembly) 

rttamAND  (  GEOMETRIC-FEASIBIUTYidecomposition,  assembly), 
STABIUJYidecomposition), 

MECHANICAL-FEASlBIUTYljiecomposition,  assembly)  ) 

end  procedure 


Figure  11:  Ptocedait  FEASIBIUTY-TEST 

5.  The  Algorithm  for  Generating  All  Assembly  Sequences 

As  discussed  in  the  previous  section,  this  research  takes  a  decomposition  approach  to  the  problem  of  generating 
assembly  sequences.  The  basic  idea  underlying  the  approach  is  to  enumerate  the  decompositions  of  the  assembly 
and  to  select  those  decompositions  that  are  feasible.  The  decompositions  are  enumerated  by  enumerating  the 
cut-sets  of  the  assembly’s  graph  of  connections.  Knowledge  of  the  feasible  decompositions  allows  the  constiuction 
of  the  ANtVOR  graph  rqnesentation  of  assembly  plans.  Each  feasible  decomposition  corresponds  to  a  hyperarc  in 
the  ANtVOR  graph  connecting  the  node  corresponding  to  the  assembly  to  the  two  nodes  corresponding  to  the  two 
subassemblies.  The  same  process  is  repeated  for  the  subassemblies  and  subsubassemblies  until  only  single  parts  are 
left 

It  has  been  shown  [12. 21]  that  the  set  of  all  cut-sets  of  a  graph  (V,E)  is  a  subspace  of  of  the  vector  space  over 
the  Galois  field  modulo  2  associated  with  the  graph.  The  vectors  in  this  vector  q»ce  are  the  elements  of  n  (£),  the 
set  of  all  subsets  of  E.  It  has  also  been  shown  that  the  fundamental  system  of  cut-sets  relative  to  a  spanning  tree  is  a 
basis  of  the  cut-set  subspace.  Therefore,  the  cut-sets  of  a  graph  can  be  enumerated  by  constructing  a  spanning  tree 
of  the  graph,  finding  the  fimdamental  system  of  cut-sm  relative  to  that  spanning  tree,  and  computing  all  the 
combinations  of  fundamental  cut-sets.  In  our  current  implemenialimi.  the  cut-sets  are  enumerated  using  a  more 
efficient  approach.  We  first  look  at  all  connected  subgraphs  having  the  cardinality  of  their  set  of  nodes  smaUer  than 
or  equal  to  half  of  the  cardinality  of  the  set  of  nodes  in  the  whole  graph.  For  each  of  these  subgraphs,  the  set  of 
edges  of  the  whole  graph  that  have  only  one  end  in  the  subgraph  defines  a  cut-set  if  their  removal  leaves  the  whole 
graph  with  exactly  two  components. 

Figure  12  shows  the  procedure  GET-FEASlBLE-DECOMPOSmONS  which  takes  as  input  the  relational  model  of 
an  assembly  and  returns  all  feasible  decompositioos  of  that  assembly.  The  procedure  first  generates  the  graph  of 
connections  fw  the  input  assembly  and  computes  the  cut-sets  of  this  graph.  Each  cut-set  corresponds  to  a 
decomposition.  The  procedure  GET-DECOMPOSmONS  is  used  to  find  the  decomposition  that  corresponds  to  a 
cut-set,  and  the  procedure  FEASIBILITY-TEST  discussed  in  the  previous  section  is  used  to  check  whether  that 
decomposition  is  feasible  or  not  The  feasible  decompositions  are  stored  in  the  list  feasible-decompositions  which 
was  empty  at  the  beginning.  After  all  cut-sets  have  been  processed,  the  procedure  returns  the  list 
feasible-decompositions. 

An  example  will  illustrate  the  computation  of  the  feasible  decompositions  of  an  assembly.  When  passed  the 
rdational  model  of  the  assembly  in  figure  3,  procedure  GET-FEASIBLE-DECOMPOSmONS  will  compute  the 
graph  of  connections  shown  in  figure  5,  and  all  its  cut-sets,  which  are  indicated  in  figure  13.  The  analysis  of  those 
cut-sets  will  indicate  the  feasible  decompositions.  The  first  cut-set  yields  a  feasible  decomposition  since  it  is 
feasible  to  join  the  cap  and  the  subassembly  made  of  the  three  other  parts.  The  second  cut-set  also  yields  a 
feasible  decomposition  because  it  is  feasible  to  join  the  subassembly  consisting  of  the  cap  plus  the  receptacle,  and 
the  subassembly  consisting  of  the  stick  plus  the  handle.  The  third  cut-set,  however,  does  not  yield  a  feasible 
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procedure  GET-FEASIBLE-DECOMPOSrnONS(assembfy) 
feasible-decomposidons  «-  NIL 
graph  GET-GRAPH-OF-CONNECnONSiassembly) 
cutsets  <-  GET-CUT SETSigraph) 
while  cutsets  is  not  empty  do 

begin  loopl 

next-cutset  <-  FIRST(cut-sets) 
cutsets  *-  TAIL(cutsets) 

next-decomposition  GET-DECOMPOSmON(next-cutset) 
if  FEASIBILITY- fESr(next-decomposition) 

then  feasible-decomposidons  UNION(feasible-decompositions,LlST{next-decomposidon)) 

end  loopl 

return  feasible-decompositions 
end  procedure 

Figure  12:  Procedure  GET-FEASIBLE-DECOMPOSITIONS 


cut-1  cut-2  cut-3 


Figure  13:  The  cut-sets  of  the  graph  of  connections  for  the  assembly  shown  in  Figure  3 


decomposition,  since  it  is  not  possible  to  join  the  stick  and  the  subassembly  made  up  of  the  three  other  parts. 
Similarly,  the  fourth  and  the  sixth  cut-sets  yield  feasible  decompositions  while  the  fifth  cut-set  does  not  Thoefore, 
procedure  GET-FEASIBLE-DECOMPOSmONS  will  return  a  list  containing  the  four  decompositions  that 
correspond  to  the  first,  second,  fourth,  and  sixth  cut-sets. 

Figure  14  shows  the  procedure  GENERATE-AND-OR-GRAPH  which  takes  the  relational  model  of  an  assembly, 
and  returns  the  AND/OR  graph  representation  of  all  assembly  sequences  for  that  assembly.  The  nodes  in  the  AND/OR 
graph  returned  are  pointers  to  relational  models  of  assemblies. 

Procedure  GENERATE-AND-OR-GRAPH  uses  the  lists  closed  and  open  to  store  the  pointers  to  the  relational 
models  of  the  subassemblies  whose  decompositions  into  smaller  subassemblies  respectively  have  and  have  not  been 
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procedure  GENERATE-AND~OR-GRAPH{assembly) 

open  «-  USTiGET-POINTmS(UST(assembly))) 

closed  <—  NIL 

hyperarcs  *-  NIL 

while  open  is  not  empty  do 

.begin  loopl 

next-subassembly  <-  FIRST(open) 
open  *-  TAIL(open) 

closed  «—  UNION(jclosed,  LIST(next-subassembly)  ) 

decompositions-of-next-subassembly  <-  GET-FEASIBLE-DECOMPOSniONS(next-subassembly) 
whUe  decomposidons-of-next-subassembly  is  not  empty  do 

begin  loop2 

next-decomposition  4-  FIRST(decompositions-of-next-subassembly) 
decomposidons-of-next-subassembly  4-  TAIL(decompositions-of-next-subassembly) 
subassemblies  4-  GET-POINTERS(next-decomposUion) 
hyperarcs  4-  UNION{hyperarcs,UST(LIST(next-subassembly,subassemblies))) 
while  subassemblies  is  not  empty  do 

begin  loop3 

next-subassembly  4-  FIRST{subassemblies) 
subassemblies  4-  TAIL{subassemblies) 

if  next-subassembly  is  not  in  open  or  in  closed,  add  it  to  open’,  otherwise  ignore  it 
end  loop3 

end  loopl 

end  loopl 

return  USTlfilosed,  hyperarcs) 
end  procedure 

Figure  14:  Procedure  GENERATE-AND-OR-GRAPH 


generated. 

The  im)cedi]re  takes  one  element  of  open  at  a  time,  moves  it  to  closed,  and  uses  procedure 
GET-FEASlBLE-DECOMPOSmONS  to  generate  all  decompositions  of  the  relational  model  pointed  by  that 
element  For  each  decomposition,  procedure  GENERATE-AND-OR-GRAPH  uses  the  procedure  GET-POINTERS  to 
get  the  pointers  to  the  relational  models  of  the  subassemblies.  Procedure  GET-POINTERS  checks  whether  each 
resulting  subassembly  has  appeared  before  or  not  If  the  subassembly  has  appeared  before,  its  pointer  is  used, 
otherwise  a  new  pointer  is  created.  The  new  pointers  are  inserted  into  open.  Each  decomposition  yields  one 
hyperatc  of  the  AND/OR  graph. 

Hgure  IS  shows  the  resulting  AND/OR  gn^h  for  the  prodiKt  shown  in  figure  3. 

A  more  efficum  implementation  of  the  method  for  the  generation  of  assembly  sequences  presented  above  will 
include  additional  tests  aimed  at  avoiding  unnecessary  computation!  One  such  test  is  to  check  whether  the 
feasibility  of  a  decomposition  follows  from  the  feasibUity  of  other  deconpositions.  For  example,  the  feasibility  of 
the  decomposition  corresponding  to  hyperaic  10  in  figure  IS  follows  firom  the  feasibility  of  the  decompositions 
corresponding  to  hyperarcs  4  and  S.  If  it  was  geometrically  and  mechanically  feasible  to  disassemble  the  handle 
from  the  whole  assembly  (hyperarc  4).  that  it  is  geometrically  and  mechanically  feasible  to  disassemble  the  handle 
fiom  a  subassembly.  And  since  the  subassembly  made  up  of  the  stick  and  the  recq)tacle  is  stable  (hyperarc  S),  it 
follows  that  the  decomposition  coneqxmding  to  hypeiarc  10  is  feasi1}le.  This  test  indicates  that  if  the 
decompositions  corresponding  to  hyperarcs  4  and  S  have  already  been  analysed  and  found  to  be  feasible,  then  it  is 
not  necessary  to  perform  the  computation  corresponding  to  procedure  FEASIBILITY-TEST  in  the  analysis  of  the 
decomposition  that  corresponds  to  hyperarc  10.  Similarly,  another  additional  test  would  check  whether  the 
unfeasibility  of  a  decomposition  follows  fix>m  the  unfeasibility  of  other  decompositions  already  analysed. 

6.  Analysis  of  the  Algorithm 

This  section  presents  an  analysis  of  the  algorithm  fot  the  generation  of  all  assembly  sequences.  First,  a  proof  of 
the  conecmess  and  completeness  of  the  algorithm  GET-FEASIBLE-DECOMPOSmONS  is  presented.  These  results 
are  then  used  to  prove  the  correctness  and  completeness  of  the  algorithm  GENERATE-AND-OR-GRAPH.  At  the 
end,  an  assessment  of  the  computation  involved  in  executing  GENERATE-AND-OR-GRAPH  is  presented. 


6.1.  The  correctness  of  algorithm  GET-FEASIBLE-DECOMPOSITIONS 

The  partial  correcmess  of  the  algorithm  GET-FEASIBLE-DECOMPOSmONS  is  immediate.  The  list  cuts  is 
initially  empty.  Only  feasible  decompositions  are  added  to  the  list  cuts.  Therefore,  the  list  returned  by 
GET-FEASIBLE-DECOMPOSITIONS  does  not  contain  any  element  that  is  not  a  feasible  decomposition  of  the 
assembly  input 

The  total  correctness  follows  from  the  fact  that  there  is  only  a  finite  number  of  cut-sets  in  a  graph.  The  list 
cut-sets  contains  irutially  all  cut-sets  of  the  graph  of  functional  cormection  corresponding  to  the  assembly  input  At 
each  execution  of  loopl,  one  elemem  is  removed  fiom  the  list  cutsets.  Therefore,  after  a  finite  number  of 
executions  of  loopl  the  list  cutsets  becomes  empty,  and  the  algorithm  terminates. 

This  proof  assumes  that  the  algcvithm  for  generating  the  cut-sets  of  a  graph  is  correct  and  complete.  As  discussed 


*Our  anrent  implaiienutiaii  consists  of  the  basic  algorithms  presented  in  the  text  and  does  not  yet  include  these  additional  tests. 
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Figure  15:  The  AND/OR  graph  for  the  assembly  shown  in  figure  3 


in  the  previous  section,  the  enumeration  of  the  cut-sets  of  a  gnq;)h  is  studied  in  graph  theory;  for  example,  Deo  [12] 
and  Liu  [21]  discuss  that  isoblem. 

This  proof  also  assumes  that  it  is  possible  to  decide  correctly  whether  a  decomposition  is  feasible  or  not,  based  on 
geometrical  and  physical  aiteiia,  as  discussed  in  the  section  4. 


62.  The  completeness  of  algorithm  GET-FEASIBLE-DECOMPOSITIONS 
There  is  a  one-to-rme  cone^tondence  between  cut-sets  in  the  gr^h  of  connections  of  an  assembly,  and  the 
decompositions  of  that  assembly.  Therefore,  since  algorithm  GET-FEASIBLE-DECOMPOSmONS  goes  over  all 
cut-sets  of  the  graph  of  connections,  all  feasible  decompositions  will  be  generated. 

As  in  the  proof  of  correctness  above,  this  proof  of  completeness  assumes  the  use  of  a  correct  and  complete 
algorithm  fm:  the  generation  of  all  cut-sets  of  a  graph,  and  a  correct  algorithm  for  deciding  the  feasibility  of  a 
decomposition. 


6  J.  The  correctness  of  algorithm  GENERATE-AND-OR-GRAPH 

List  closed  is  updated  at  only  one  point,  and  it  only  gets  dements  that  were  previously  in  the  open  list  The  open 
list  contains  initially  a  pointer  to  the  relational  modd  of  the  assembly  input,  which  is  a  node  of  the  AND/OR  graph. 
List  open  is  updated  inside  loop3  where  it  gets  pointers  to  the  relational  models  of  the  subassemblies  that  are  pan  of 
a  feasible  decomposition,  and  therefore,  are  nodes  of  the  ANO/OR  grq>h.  Thnefore,  the  elements  in  the  open  list, 
and  consequently  the  dements  in  the  closed  list,  are  always  pointers  to  relational  models  either  of  the  original 
assembly,  or  of  subassemblies  that  take  part  of  a  feasiUe  decompositiem. 

The  hyperarcs  hst  is  initially  empty.  It  is  updated  only  inside  loop2  where  it  gets  the  hyperarc  corresponding  to  a 
feasible  decomposition.  Therefore,  algorithm  GET-FEASIBLE-DECOMPOSmONS  can  only  return  a  set  of  nodes 
and  a  set  of  hyperarcs  of  the  AND/OR  gnph.  This  establishes  the  partial  correctness  of  the  algorithm. 

List  open  gets  only  subassemblies  and  no  subassembly  is  inserted  more  than  once.  Since  thoe  is  a  finite  number 
of  subassemblies,  the  algorithm  terminates.  This  establishes  the  total  cotreemess  of  the  algorithm. 


6.4.  The  completeness  of  algorithm  GENERATE-AND-OR-GRAPH 
Since  algorithm  GET -FEASIBLE-DECOMPOSITIONS  is  complete,  all  possible  decompositions  of  all 
subassemblies  that  are  inserted  into  the  list  open  yield  a  hyperaic.  Furthermore,  all  subassemblies  that  result  from  a 
decomposition  are  inserted  into  list  open,  and  later  are  moved  to  list  dosed.  Therefore,  the  first  list  returned 
contains  all  subassemblies  that  resulted  from  some  decomposition,  and  the  second  list  returned  contains  one 
hyperarc  for  each  decomposition  of  each  subassembly. 


6.5.  Complexity 

The  amount  of  computation  involved  in  the  generation  of  the  AND/OR  griq)h  for  a  given  assembly  depends  on  the 
number  N  of  parts  that  make  iq)  the  assembly,  on  how  interconnected  those  parts  are,  and  also  on  the  resulting 
AND/OR  graph. 

The  number  of  prospective  decompositions  (i.e.  cut-sets  of  the  graphs  of  functional  connections)  that  must  be 
analysed  will  be  used  in  this  section  as  a  measure  of  the  ^ount  of  computation  involved  in  the  generation  of  all 
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assembly  sequences^.  Two  models  for  how  the  parts  in  the  assembly  are  interconnected  are  considered  in  order  to 
im)vide  bounds  in  the  estimate  of  computational  complexity: 

1.  a  strongly  connected  assembly  in  which  every  part  is  connected  to  every  otho'  part;  and 

2.  a  weakly  connected  assembly  in  which  there  are  A/-1  connections  between  the  N  parts,  with  the 
connection  being  between  part  the  and  the  ({>!)'’*  paits. 

And  three  possibilities  for  the  resulting  AND/OR  graph  are  c(Hisidere± 

1.  a  balanced  tree  AND/OR  graph  in  which  there  is  at  most  one  hyperarc  leaving  each  node  and  this 
hypeiarc  points  to  two  nodes  whose  coiesponding  subassemblies  either  have  the  same  number  of 
parts,  or  their  number  of  parts  differ  by  one; 

2.  one-part-at-a-dme  tree  AND-OR  graph  in  which  there  is  at  most  one  hyperarc  leaving  each  node,  and 
•this  hypetarc  points  to  two  nodes  one  of  which  coneqxmds  to  a  one-part  subassembly;  and 

3.  a  network  AND/OR  grtq)h  in  which  there  are  as  many  hypearcs  leaving  each  node  as  there  are  cut-sets  in 
the  graph  of  functional  crauiections  al  the  node’s  corresponding  subassembly. 


The  resulting  total  number  D  of  decompositions  that  must  be  analysed  as  a  function  of  the  number  N  of  parts  that 
make  up  the  assembly  for  each  possible  combination  of  how  the  parts  are  interconnected  and  the  type  of  the 
resulting  AND/OR  graph  is: 

1.  Weakly  coimected  assemblies: 

a.  Balanced  tree  AND/OR  grs^h:  the  number  of  prospective  decompositions  that  must  be  analysed 
is  W-1  fcv  the  initial  assembly,  N-2  for  all  subassemblies,  N-4  for  all  subsubassemblies,  and  so 
on.  Therefore^, 


f)»(A/-l)+(A/-2)+(iV--4)+  •••  +(Ar-2 


iiit(log2  A/) 


)  = 


.  [im(log2N)+I] 

=  y  (JV-2)  =  A/-[int(log2A/)  +  l]-2  ^  +1 

b.  One-part-at-a-time  tree  AND/OR  grq}h:  the  luimber  of  proq)ective  dectunpositions  that  must  be 
analysed  is  N-l  for  the  initial  assembly,  N~2  for  the  (^-2)  part  subassembly,  N-3  for  the 
(A/-2)-part  subassembly,  and  so  on.  Therefrae, 

£)-(A/-l)+(A/-2)+(/V-3)+  •••  +2+1  «  (^-0  » 

c.  Network  AND/OR  graph:  the  number  of  prospective  decompositions  that  must  be  analysed  is 
N-l  for  the  A/-part  subassembly,  N-2  for  each  of  the  two  (/V-l)-part  subassemblies,  N-3  for 
each  of  the  three  (A/-2}-part  subassemblies,  and  so  on.  Therefore, 

£)=  l(A/-l)  +  2(N-2)  +  3(A/-3)+  •••  +(A/-1)1  = 


N-l 


iN+l)-N-iN-l) 

6 


ovenU  complexity  of  algorithm  CENERATE-AND-OR-GRAPH  ihoold  lake  iitto  account  the  computation  invcdved  in  generating  the 
cut-sets  of  the  graph  of  functional  connections. 

^e  use  the  notation  int(x  }  to  represent  the  largest  integer  that  it  less  than  or  equal  to x.  For  example,  int  (3  )= 3  and  int  ( 3.5  )  =  3. 
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2.  Strongly  connected  assemblies; 

a.  Balanced  tree  AND/OR  grig)h:  the  number  of  prospective  decompositions  that  must  be  analysed 

is  (2  -1)  for  the  initial  assembly,  (2  ^  -1)  +  (2  *  - 1 )  for  all  subassemblies, 

toO  in.a 

(2  ^  -1)  +  (2  -1)  +  (2  <  -1)  +  (2  “  -1)  for  all 


subsubassemblies,  and  so  on.  Therefore, 
imflogjA/)  2i-i 

b.  One'part-at-a-time  tree  AND/OR  grt^:  the  number  of  prospective  decompositions  that  must  be 
N-l  N-2 

analysed  is  (2  - 1)  for  the  ^-part  subassembly,  (2  - 1 )  for  the  (iV- 1  )-patt  subassembly, 

iV-3 

(2  -1)  for  the  (iV- 2  >part  subassembly,  andso  on.  Therefore, 

N-l  N-2  N-i  N 

D=(2  -l)+(2  -l)+(2  -1)+ ••• +(2-l)=  2  -iV-l 


c.  Network  AND/OR  graph:  the  number  of  proq)ective  decompositions  that  must  be  analysed  is 

N-l  N-2  ^  «  V 

(2  -1)  for  the  N-part  subassembly,  (2  -1)  for  each  of  the  (^i)  CN-l)-part 

N-2  ,  „  ^ 

subassemblitt,  (2  -1)  each  of  the  (^_2)  (N--2)-part  subassonblies,  and  so  on. 

Therefore, 


For  each  of  the  three  possibilities  (rf  the  resulting  AND/OR  graph,  table  2  shows  the  number  of  decompositions  that 
must  be  analysed  for  weakly  connected  assemblies  and  table  3  shows  the  number  of  decnnpositions  that  must  be 
analysed  for  strongly  connected  assemblies,  as  a  function  of  the  number  of  parts  that  make  op  the  product  The 
figures  in  table  3  ate  given  as  a  reference  since  it  is  very  unlikely  that  there  would  be  a  twenty-part  assembly  in 
which  every  part  is  connected  to  every  other  part 

The  results  above  take  into  account  the  fact  that  the  type  of  the  resulting  AND/OR  grq)h  is  not  known  a  priori.  For 
example,  for  the  weakly  connected  assembly  whose  AND/OR  graph  is  a  balanced  tree,  all  the  N-l  cut-sets  of  the 
whole  assembly  were  included  in  the  number  of  decompositions  that  are  tested,  although  thoe  is  only  one  cut-set 
that  yields  two  subassemblies  that  have  the  same  number  of  parts. 

As  discussed  in  the  end  of  section  S,  a  mote  efficient  implementation  of  the  method  for  the  generation  of 
assembly  sequences  presented  in  this  paper  will  include  additional  tests  aimed  at  avoiding  unnecessary  computation. 
One  such  test  is  to  check  whether  the  feasibility  of  a  decomposition  follows  from  the  feasibility  of  other 
decompositions.  In  the  case  of  strongly  connected  assemblies  in  which  all  decompositions  of  all  subassemblies  are 
feasible,  the  computation  can  be  significantly  reduced  if  this  test  is  performed  before  analysing  each  decomposition. 
Since  all  decompositions  of  the  whole  assembly  are  feasible,  all  decompositions  of  all  subassemblies  should  also  be 
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Table  2:  The  number  of  decompositions  that  must  be  analysed  for  each  type  of 
resulting  AND/OR  gn^h,  as  a  function  of  the  number  of  parts,  for 
weakly  connected  assemblies. 


Number  of  Parts 

N 

Balanced-tree 
AND/OR  graph 

One-part-at-a-time 
AND/OR  grtqih 

network 

AND/OR  graph 

2 

1 

1 

1 

3 

3 

3 

4 

4 

5 

6 

10 

5 

8 

10 

20 

6 

11 

15 

35 

7 

14 

21 

56 

8 

17 

28 

84 

9 

21 

36 

120 

10 

25 

45 

165 

15 

45 

105 

560 

20 

69 

190 

U30 

25 

94 

300 

2,600 

30 

119 

435 

4.495 

Table  3:  The  number  of  decompositions  that  must  be  analysed  for  each  type  of 
resulting  AND/OR  graph,  as  a  fiinctiai  of  the  number  of  parts,  for 
strongly  connected  assemblies. 


Number  of  Parts 

N 

Balanced-tree 
AND/OR  graph 

One-part-at-a-time 
AND/OR  graph 

network 

AND/OR  graph 

2 

1 

1 

1 

3 

4 

4 

6 

4 

9 

11 

25 

5 

20 

26 

90 

6 

39 

57 

301 

7 

76 

120 

966 

8 

145 

247 

3,025 

9 

284 

502 

9330 

10 

551 

1,013 

28.501 

15 

16,604 

32,752 

7.141,686 

20 

525389 

1.048355 

1.742343,625 

25 

16,783,550 

33354,406 

423,610,750390 

30 

536304.119 

1,073,741,793 

102,944,492,305301 
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feasible.  Therefoie,  a  suable  additicmai  test,  the  total  number  of  decompositions  that  must  be  analysed  is 
reduced fiom -  2  to(2  -1). 


7.  Conclusion 

A  ccHiect  and  complete  algorithm  for  the  generation  of  all  mechanical  assembly  sequences  was  presented.  To  the 
authtxs*  knowledge,  no  inevioos  algorithm  for  the  generation  of  all  mechanical  assembly  sequences  has  been  proven 
correct  and  complete. 

The  problem  of  generating  assembly  sequences  was  transformed  into  the  equivalent  problem  of  generating 
disassembly  sequences.  The  algorithm  operation  consists  in  looking  at  all  the  decompositions  of  the  assembly,  i.e. 
all  the  ways  the  assembly  can  be  split  into  two  subassemblies.  This  is  done  by  genoating  all  cut-sets  of  the 
assembly’s  graph  of  connections,  and  checking  which  cut-sets  correspond  to  feasible  decompositions.  A 
decomposition  is  feasible  if  it  possible  to  obtain  the  assembly  by  jtMning  the  two  subassemblies.  The  same  process 
is  repeated  for  the  subassemblies,  for  the  subsubassemblies,  and  so  on,  until  only  single  parts  are  left.  At  the  end, 
the  AND/OR  graph  rqnesentaticMi  of  assembly  sequences  is  returned. 

The  algorithm  also  lends  itself  to  an  interactive  implementadcm  in  which  a  computer  program  generates  questions 
that  are  answered  by  a  human  expert  Each  quesdai  addresses  the  feasibility  of  a  decomposition.  But  unlike 
previous  methods  [6. 11],  it  is  possible  to  have  a  computer  irogram,  instead  of  a  human,  to  answer  the  questions 
directly  from  a  descrqKion  d  the  assembly.  Our  current  imptementation,  which  has  the  restrictions  on  the  types  of 
assemblies  discussed  in  section  3,  includes  {Hograms  that  answv  the  questions. 

An  approach  to  compute  the  answer  to  the  question  of  whether  it  is  feasible  to  obtain  a  given  assembly  by  joining 
two  subassemblies  was  (nesented.  This  approach  is  based  on  the  use  of  a  relational  model  description  of  the 
assembly.  The  model  includes  three  types  of  entities:  parts,  contacts,  and  attachments;  it  also  includes  a  set  of 
relationships  between  entities.  Both  entities  and  relationships  can  have  attributes.  To  decide  whether  a  given 
decomposition  is  feasible,  three  predicates  must  be  computed,  using  the  data  in  the  relational  model: 

•  The  GEOMETRIC-FEASIBIIUY  predicate  which  is  true  if  there  exists  a  collision-ftee  path  to  bring  the 
two  physical  subassemblies  into  contact  fiom  a  situation  in  which  they  are  sufficiently  far  spirt 

•  The  MECHANICAL-FEASIBILITY  predicate  which  is  true  if  it  is  feasible  to  establish  the  attachments 
that  act  on  the  contacts  of  the  decomposition. 

•  The  STABILITY  predicate  which  is  true  if  the  parts  in  each  subassembly  maintain  their  relative  position 
and  do  not  break  contaa  spontaneously. 

The  key  assumption  in  proving  the  conectness  of  the  algorithm  was  that  it  is  always  possible  to  decide  correctly, 
based  on  geometrical  and  physical  criteria  (i.e.  using  the  three  predicates  above),  whether  it  is  feasible  to  obtain  a 
given  assembly  by  joining  two  subassemblies. 

The  amount  of  computation  involved  in  generating  aU  mechanical  assembly  sequences  was  assessed  by 
determining  the  number  of  decompositions  that  must  be  analysed.  That  amount  dq)ends  not  only  on  the  number  of 
parts  and  on  how  they  are  interconnected,  but  <m  die  solution  AND/OR  gnqih  as  well.  The  least  amount  of 
computation  occurs  for  weakly  connected  assonblies  in  which  each  subassembly  has  only  one  feasible 
decomposition  and  that  decomposition  yields  two  subassemblies  whose  number  of  parts  are  either  equal  or  differ  by 
one.  The  maximum  amount  of  computation  occurs  for  strongly  connected  assemblies  in  which  all  decompositions 
of  all  subassemblies  are  feasible.  This  worst  case,  howevo*,  is  very  unlikely  to  occur  in  practice.  Furthermore, 
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additional  simple  tests  discussed  in  section  S  can  reduce  the  amount  of  computation. 


In  practice,  an  evaluation  of  the  alternative  assembly  sequences  generated  by  the  algorithm  presented  in  this  paper 
is  needed  in  order  to  choose  the  sequence  that  will  be  tK;tually  used  in  the  assembly  process.  Different  evaluation 
fiintions  have  been  explored  including  a  function  based  on  parts  entropy  [32, 33],  and  a  function  based  on  the 
complexity  of  assembly  tasks  and  the  stability  of  subassemblies  [17]. 

It  is  also  possible  to  implement  an  interactive  system  in  which  a  computer  program  generates  the  alternative 
sequences,  as  described  in  this  paper,  and  a  human  expert  then  selects  the  best  one.  Still  another  possibility  would 
be  to  use  an  evaluation  function  for  a  preselection  of  "good”  alternative  sequences  and  that  have  a  human  expert  to 
make  the  final  choice. 

Whenever  the  amount  of  computation  exceeds  the  availid>Ie  computational  resources,  at  least  two  strategies  may 
be  followed: 

1.  The  number  of  parts  can  be  artificially  reduced  by  treating  subassemblies  as  single  parts.  An  analysis 
of  the  graph  of  connections  may  indicate  the  clusterings  of  parts  that  yield  bigger  reductions  in  the 
amount  of  computation. 

2.  The  algorithm  generates  fewer,  hopefully  the  best,  sequences  using  some  heuristics  to  guide  the 
generatkm  of  assembly  sequence.  Such  heuristics  should  be  compatible  with  the  evaluation  function 
used  to  choose  among  the  alternative  assembly  sequences. 

In  bodi  strategies,  the  computation  will  be  reduced  stt  the  expense  ci  the  completeness,  since  not  all  possible 
sequences  will  be  generated.  The  devolpment  of  a  procedure  to  cluster  parts  into  subassemblies  to  obtain  a 
hierarchical  model  of  the  assembly,  and  the  development  ci  good  heuristics  to  guide  the  generation  of  assembly 
sequences  are  issues  for  future  research. 
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Appendix 


L  Reasoning  about  the  Feasibility  of  Local  Translations  for  Robotic  Assembly  of  a  Part 
Constrained  by  Planar  Contacts 

LI.  Introduction 

The  high  level  of  planning  the  assembly  of  a  {soduct  can  be  viewed  as  a  path  search  in  the  state-space  of  possible 
configurations  of  the  set  of  parts  that  comprises  the  product  [17].  The  initial  state  corresponds  to  the  configuration 
in  which  all  parts  are  discoimected  from  each  other.  The  goal  state  corresponds  to  the  configuration  in  which  the 
parts  ate  properly  joined.  The  moves  correspond  to  the  assembly  operations,  since  they  change  one  state  into 
another. 

A  rpinpiftte  description  of  die  product  is  available  for  planning  purposes.  This  description  includes  the  shape  of 
the  parts,  thdr  relative  positions,  and  the  spatial  and  mechanical  relations  between  parts. 

The  search  can  be  conducted  backwards  fiom  the  goal  state  to  the  initial  state.  The  moves  in  the  backward  search 
correspond  to  disassembly  tasks  which  are  defined  to  be  the  reverse  of  fisasible  assembly  tasks.  The  preconditions 
fw  a  disassembly  task  [34]  include: 

1.  release  of  attachments. 

2.  stability  of  subassemblies. 

3.  separability  of  subassemblies: 

a.  local  analysis  •  test  incremental  motion; 

b.  global  analysis  •  find  global  trajectory. 

The  local  analysis  consists  of  checking  whether  there  exists  an  incremental  motion  of  one  part  or  subassembly 
that  is  not  blocked  by  any  one  of  its  contacts  with  other  parts.  For  many  types  of  contacts  there  are  very  few  feasible 
motions  between  the  parts.  Fot  example,  a  cylindrical  pin  in  a  hole  can  either  translate  in  the  direction  of  the  axis,  or 
rotate  around  the  axis.  Whenever  the  part  or  subassembly  under  consideration  has  such  a  constrairung  contact,  the 
local  analysis  can  be  performed  by  checking  whether  at  least  one  of  the  few  motions  that  are  compatible  with  the 
most  restrictive  contact  is  also  compatible  with  all  other  contacts. 

The  local  analysis  is  mcne  difficult  when  the  part  (<^  subassembly)  to  be  disassembled  is  constrained  by  planar 
contacts  only.  Each  planar  contact  leaves  an  infinite  number  of  unconstrained  directions  along  which  translation  is 
possible.  In  this  case,  the  intersection  of  the  sets  of  translations  that  are  not  blocked  by  each  contact  caruiot  be  found 
by  discrete  search  over  a  finite  set  of  directions. 

This  appendix  presents  an  efficient  procedure  to  obtain  explicitly  the  set  of  directions  along  which  an  object  that 
is  constrained  by  several  planar  contacts  can  translate.  In  addition  to  answering  the  question  of  whether  there  is  a 
direction  along  which  translation  is  feasible,  the  procedure  also  produces  a  representation  of  the  set  of  all  those 
directions. 
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12.  Background 


Shoham  [35]  and  Mani  [26]  analyzed  the  freedom  of  two  dimensional  objects  when  in  contact  with  other  objects. 
Although  both  Studies  include  translational  and  rotational  freedom,  they  are  restricted  to  two  dimensional  objects, 
and  neither  one  indicates  how  the  analysis  can  be  scaled  up  to  three  dimensions. 

lain  and  Donath  [19]  analyzed  the  translational  and  rotatkmal  freedom  of  parts  in  an  assembly  but  with  the 
constraint  that  the  existing  contacts  between  parts  cannot  be  broken.  They  do  not  show  how  their  aqpproach  can  be 
extended  to  deal  with  the  case  of  breaking  of  contacts. 


Ejiri  et  al.  [14]  pttqiosed  the  use  of  restraint  vectors  to  decide  whether  a  part  constrained  by  planar  contacts  could 
translate.  The  restraint  vecttv  of  a  part /’was  defined  to  be 

®+l* 

where 

^  _  r  1  if  /*  is  restrained  along  direction  r 
''~\0  if  /*  is  not  restrained  along  direction  r 


and  the  logical  possibility  of  disassembling  part  P  was  decided  by  the  logical  formula 


a  €L 
“+jr  -jc 


_  r  1,  when  possible 
~\0,  when  impo^ble. 


(1) 

That  formula  corresponds  to  the  r^dtement  thru  the  part  be  free  in  the  upper  (positive  z)  direction  and  both  x  and  y 
are  free  in  either  the  positive  or  the  negative  direction.  Although  the  above  togical  formula  is  a  sufficient  condition, 
it  is  not  necessary.  It  is  also  not  difficult  to  think  of  a  situation  in  which  disassembling  a  part  is  feasible  and  the 
logical  formula  yields  0.  Hgtne  16  shows  one  such  situation. 


ap=(1, 1,0, 1,0,1) 


Figure  16:  Part  P  can  move  but  the  logical  formula  (1)  yields  0 
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Within  the  wtnic  in  Idnematics,  Asada  and  By  [3]  intnxhiced  the  concept  of  Automatically  Recor^gured  Fixturing 
which  is  a  fixturing  system  that  can  be  adig)ted  to  hold  difierent  workpaits.  It  consists  of  a  number  of  fixture 
elements  that  can  be  placed  on  a  flat  hcmzontal  table  to  conform  to  the  geometry  of  the  woilcpan.  to  be  fixtured.  The 
table  has  magnetic  chucking  capability  which  can  be  activated  to  secure  the  fixturing  elements  in  place.  When 
completely  fixtured,  the  w(xlq;)art  will  be  in  contact  with  a  number  of  fixturing  elements  that  will  constrain  its 
movements  completely.  The  {socess  of  fixturing  one  woikpart  starts  with  a  positioning  phase  in  which  the  wcxkpait 
is  brought  into  contact  with  a  subset  of  the  fixturing  elements  which,  in  this  p^)er,  will  be  referred  to  as  guiding 
elements.  Once  positioning  is  achieved,  addititmal  elements  are  placed  on  the  table  to  constrain  the  woriq>art 
completely. 

Asada  and  By  carried  out  a  kinematic  analysis  aimed  at  answering  the  following  questions: 

1.  Is  the  location  of  the  woriqnrt  that  achieves  contact  with  all  guiding  elements  unique? 

2.  Is  this  location  that  achieves  contact  with  the  guiding  elements  accessible/detachable? 

3.  Do  the  additional  elements  (togedier  with  the  guiding  elements)  constrain  the  worirpait  completely? 

The  last  two  questions  are  similar  to  the  problem  addressed  in  this  paper.  In  the  second  question,  they  are 
interested  in  guaranteeing  that  there  are  feasible  local  motions  for  the  workpart,  so  it  can  be  brought  into  contact 
with  all  guiding  elements.  In  the  third  question,  they  are  interested  in  guaranteeing  that  there  are  no  feasible  local 
motions  (Le.  that  the  woriqmrt  is  constrained  completely). 

Asada  and  By  modeled  die  contacts  between  the  wtxkpart  and  the  fixturing  elements  as  point  contacts,  and 
derived  conditions  fcv  the  feasibility  of  local  motions,  including  both  translations  and  rotations.  Those  constraints 
were  used  to  check  whether  the  configuration  of  the  fixturing  elements  would  constrain  the  workpart  completely. 
They  did  not  address  how  to  determine  the  set  of  incremental  motions  that  satisfy  the  derived  conditions. 

More  recent  research  on  robotic  planning  [22, 23]  has  aimed  at  enabling  robots  to  execute  tasks  specified  in 
task-level  commands  such  as 

move  <part-id>  to  <location-specificcttion> 

in  which  the  second  term  within  angle  brackets  specifies  a  configuration  (a  position  and  an  orientation)  either  as  a 
homogeneous  transform  matrix  (v  as  a  set  of  spatial  relationsh^  among  objects.  The  translation  of  a  task-level 
command  into  robot-level  commands  involves  selecting  fixtures,  grasping  points,  gross  motions,  fine  motions,  etc. 

It  is  clear  that  procedures  that  are  able  to  construct  a  path  for  a  part  from  an  initial  configuration  to  a  final 
configuration  can  also  be  used  to  answer  whether  there  exists  a  direction  in  which  local  translation  of  that  part  in  the 
initial  configuration  is  feasible.  If  there  is  a  path,  there  is  a  direction  in  which  local  translation  is  feasible. 
Lozano-Pdrez  [24]  shows  one  procedure  to  construct  a  path  that  avoids  obstacles  and  lists  the  most  significant 
literature  on  that  subject  The  procedures  to  construct  a  path,  however,  involve  extensive  computation,  and  therefore 
their  use  in  the  high  level  of  planning  will  weaken  the  planner  efficiency.  One  of  the  major  advantages  of 
hierarchical  plaruiing  [31]  is  the  possibility  of  abstracting  the  details  at  the  high  level. 

L3.  Representation  of  Local  Constraints 

In  most  cases  of  two  parts  or  subassemblies  in  contact,  a  pure  rotation  of  one  with  respect  to  the  other  will  not 
separate  them.  In  these  cases  the  motion  must  include  a  nonzero  translational  componoit  in  order  to  separate  parts 
in  contact  Therefore,  to  decide  whether  two  parts  in  contaa  can  separate  fiom  each  other,  the  local  analysis  can 
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focus  on  translational  motions  only^. 

Let  P  and  Q  be  two  parts  that  have  one  planar  contact  Let  n  be  a  vector  perpendicular  to  the  contact  plane,  and 
pointing  towards  part  P.  Part  Q  blocks  translations  of  part  P  that  have  negative  projections  over  n .  Therefore,  to 
decide  whether  part  P  can  translate  by  a  vector  / ,  it  is  necessary  to  check  whether  t  •  n  (the  scalar  product  of  t_ 
and  n  )  is  greater  than  at  equal  to  zero. 

In  general,  a  part  P  has  ^  planar  contacts  with  other  pans:  let  be  a  vector  perpendicular  to  the  plane  of  the  i 
planar  contact,  pointing  towards  pan  P.  Then,r  must  satisfy 

^0  1=1,2,  •••.iV  (2) 

in  order  to  be  a  feasible  translation  for  pan  P. 

The  set  of  inequalities  (2)  is  a  necessary  but  not  sufficient  condition  for  the  global  translation  t  of  pan  P,  since 
other  parts  that  are  not  directly  in  contact  with  P  may  also  constrain  its  movements.  For  local  analysis,  distant 
objects  do  not  interfere,  and  system  (2)  becomes  a  sufficient  condition.  Moreover,  if  satisfies  system  (2),  so  do 
all  vectors  yt^  for  any  scalar  y  greater  than  zero;  and  it  is  always  possible  to  pick  y  sufficiently  small  to  guarantee 
that  a  translation  by  of  pan  P  is  feasible.  Therefore,  to  answer  whether  pan  P  can  translate  locally,  it  is 
sufficient  to  answer  whether  system  (2)  has  a  nonzero  solution. 

Each  inequality  in  system  (2)  divides  the  space  into  two  halfspaces.  The  set  of  vectors  satisfying  the  system  of 
inequalities  Q),  which  is  the  intersection  of  finitely  many  halfspaces,  is  a  polyhedral  convex  cone.  Polyhedral 
convex  cones  may  have  several  different  shtq)es  and  enumeration  of  these  shapes  will  be  useful  for  the  search 
procedure.  Figure  17  shows  one  example  of  a  polyhedral  convex  cone  which  is  the  intersection  of  five  halfspaces, 
each  one  defined  by  a  plane  that  goes  through  the  origin  and  that  is  perpendicular  to  a  vector  n,- .  Therefore,  the  cone 
can  be  characterized  as  the  set  of  vectors  that  have  greater  than  or  equal  to  zero  projection  ova-  vectors 
a  j  ,  ^2  ,  /I3  ,  n  4  ,  n  j  ,  which  are  perpendicular  to  the  five  faces  of  the  cone,  and  have  the  ^ipropriate  (i.e.  towards 
the  inside)  orientations.  Alternatively,  the  same  polyhedral  convex  cone  can  be  characterized  as  the  set  of  positive 
linear  combinations  of  vectors  e  j .  £2  •  £3  •  £4  •  £5  •  which  have  directions  along  the  five  edges  of  the  cone,  and  the 
appropriate  orientations  [16]. 

These  alternative  rqnesentations  of  the  polyhedral  convex  cone  and  their  properties  may  be  defined  as  follows: 

Definition  5;  Given  a  polyhedral  convex  cone  C,  any  set  of  vectors  F=  { v, ,  V2  ,  •  ■  •  ,  Vy  )  with  the 
property  that  any  vector  x €  C  has  positive  projection  over  all  vectors  v,- e  V  (i.e.  x  •  v-^0  for 
i=  1,2,  •  •  •  ,7),  is  called  a  tangential  representation  of  C. 

Definition  6:  Given  a  polyhedral  convex  cone  C,  any  set  of  vectors  (e  i  >  £2  ,  ■  ■  -  ,  £y  ]  such  that 
any  positive  linear  combination  of  the  e  j  ,  £2  .  •  •  •  .  £y  yields  a  vector  in  C  and,  convosely,  any  vector 
X  in  C  can  be  expressed  as  a  positive  linear  combination  of  £  | ,  £2  .  ■  -  ■  ,  £y  (i.e.  x  =  a,-  £,•  with 

a,'  2  0  1=1,2,  •  •  •  ,  7,  if  and  only  if  xe  Qis  called  a  point  representation  of  C. 

Definition  7:  Two  distinct  (point  or  tangential)  representations  of  the  same  polyhedral  convex  cone  are 
said  to  be  equivalent  representations. 

Definition  8:  A  tangential  representation  of  a  polyhedral  convex  cone  is  said  to  be  a  minimal 
tangential  representation  if  it  has  no  equivalent  tangential  rqnesentation  with  fewer  vectors. 


^Theie  are  caset  where  two  pant  in  contaa  can  only  separate  from  each  other  if  one  undergoes  a  pure  rotation  followed  by  a  translation.  These 
cates  correspond  to  mom  complex  contacts,  and  therefore  require  the  use  of  mote  complex  models.  In  products  designed  for  assemby  [2,  S], 
however,  these  mote  conqslex  contacts  are  rate. 
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Figure  17:  A  polyhedral  convex  cone  which  is  the  intersection  of  five  halfspaces 


Definition  9:  A  point  representation  of  a  polyhedral  convex  cone  is  said  to  be  a  minimal  point 
representation  if  it  has  no  equivalent  point  representation  with  fewer  vectors. 

The  minimal  point  representation  of  polyhedral  convex  cones  is  used  in  the  local  analysis  of  feasible  directions  of 
motion  because  one  can  readily  check  if  the  set  of  feasible  nonzero  translations  is  empty  by  checking  whether  the 
point  rq}resentation  has  a  nonzero  vector.  It  is  also  useful  as  a  basis  for  the  global  analysis  because  it  allows  the 
enumeration  of  the  feasible  translations. 

In  the  next  section,  a  syntax  for  a  computer  representation  of  polyhedral  convex  cones  in  is  developed  based 
on  the  formulation  described  above.  A  procedure  is  defined  which  finds  the  computer  representation  firom  a 
tangential  representation.  The  procedure  provides  the  basis  for  testing  the  feasibility  of  disassembly  operations. 

L4.  Search  Procedure  for  Feasible  Local  Translations 

The  solution  procedure  is  based  on  enumeration  of  the  10  possibilities  for  the  shape  of  a  polyhedral  convex  cone 
in  the  three  dimensional  space  R^.  These  are  listed  in  table  4.  A  syntax  for  a  computer  representation  of  polyhedral 
convex  cones,  which  incorporates  explicitly  their  shapes  is  shown  in  figure  18.  This  representation  is  compact,  yet 
captures  all  the  information  needed  in  the  procedure  to  find  the  solutions  of  the  system  of  inequalities  (2). 

Figure  19  shows  procedure  SOLVE  which  takes  as  input  a  tangential  representation  of  a  polyhedral  convex  cone, 
H . «2  •  ■  ■  ■  'ILn  }  which  need  not  be  minimal,  and  returns  its  computer  representation,  in  the  syntax  shown 
in  figure  18. 
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Table  4:  The  possible  shapes  of  polyhedral  convex  cones  in  three  dimensional  space 


shape 

definition 

comments 

SPACE 

All  points  in  R^. 

This  is  a  degenerate  case  in  which  iV  in  the 
system  of  inequalities  (2)  (Le.  the  number  of 
planes)  is  zero. 

HALFSPACE 

All  the  points  on  one  side  of  a  plane  and  the 
points  on  the  plane. 

Typically,  this  is  the  case  in  which  N  in  system 
(2)  is  one,  but  it  is  also  the  case  in  which  all 
vectors  R  i  >  n  2 .  -  -  -  ,  are  parallel  and  have 

the  same  orientation. 

QUADRANT 

The  points  within  the  intersection  of  two 
halfspaces  whose  defining  planes  are  not 
parallel. 

Typically,  this  is  the  case  in  which  N  is  two,  and 
the  two  vectors  n  j  and  £2  are  not  parallel. 

POLYGONAL 

The  points  within  the  intersection  of  three  or 
nxMe  halfspaces  when  no  plane  exists  that 
contains  all  points  in  the  cone. 

This  is  the  shape  of  the  cone  in  figure  17. 

PLANE 

All  the  points  that  lie  on  a  plane  that  goes 
through  the  origin. 

This  is  the  two  dimensional  correspondent  of 
SPACE.  Typically,  this  is  the  case  in  which  N  is 
two,  and  the  two  vectors  are  parallel  and  have 
opposite  orientations. 

HALFPUNE 

All  the  points  that  lie  on  a  plane  that  goes 
through  the  origin,  and  to  one  side  of  a  line  on 
that  plane  that  also  goes  throug}i  the  origin. 

Typically  this  is  the  case  in  which  N  is  three, 
two  of  the  vectors  are  parallel  and  have  opposite 
(mentations,  and  the  third  vector  is  not  parallel 
to  the  other  two. 

SECTOR 

All  the  points  that  lie  on  a  plane  th^  goes 
through  the  origin,  and  also  to  one  side  of  two 
lines  on  that  plane  that  also  go  through  the 
origin. 

Typically,  diis  is  the  case  in  which  N  is  four, 
two  of  the  vectors  are  parallel  and  have  opposite 
orientations,  and  both  the  other  two  vectors  are 
not  parallel  to  any  other  vector. 

LINE 

All  the  points  on  one  straight  line. 

One  example  of  this  is  the  case  in  which  N  is 
three  and  the  three  vectors  lie  on  a  plane,  with 
no  two  vectors  parallel,  and  no  one  of  the  three 
vectors  can  be  expressed  as  a  positive  linear 
combination  of  the  other  two. 

HALFLINE 

All  the  points  on  one  straight  line,  to  one  side  of 
the  origin. 

Typically,  this  is  the  case  in  which  N  is  five, 
four  vectors  define  a  line,  and  a  fifth  vector  is 
not  parallel  to  any  one  of  the  other  four. 

POINT 

The  origin. 

This  is  also  a  degenerate  case  in  which  the  only 
solution  to  the  system  of  inequalities  is  the  zero 
vector. 
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<con«>  ■’  (SPACE)  I  (HALFSPACE  <v*ctor»  | 

(QUADRANT  (<'7«ctor>  <vactor>)  )  | 

(POLYGONAL  (<Biln-taxxgttnt-repr>  <mln-polnt-repr>)  )  | 
(PLANE  <‘vector>)  |  (HALFPLANE  (<vector>  <vector>)  )  | 
(SECTOR  <mln-polnt-rapr>)  |  (LINE  <vector>)  | 

(HALFLINE  <v«ctor>)  |  (POINT) 

<mln-tangant-repr>  «  (<v«ctor-sequenca>) 

<inln-polnt-r«pr>  a  (<v«ctor-saqu«nea>) 

<v«ctor-sequ«nc«>  «  <vactor>  |  <vactor>  <vector-sequ«nca> 


Figure  18:  The  computer  representations  of  cones 


procedure  SOLVE(tanrep) 
solution  *-  (SPACE) 

while  FIRST(solution)  ^  POINT  and  teuirep  is  not  empty  do 


begin 

n  <-  FIRST(tanrep) 
tanrep  <-  TAIL(tanrep) 
solution  <-  INTER(solution  n  ) 

end 

return  solution 
md  SOLVE 


Figure  19;  The  procedure  SOLVE 


The  solution  procedure  consists  of  finding,  successively,  a  computer  representation  for  the  sequence  of  cones 
Co.CpCj,  •  •  •  ,Cn-  The  polyhedral  convex  cone  Cq  is  the  whole  space  R^,  and  the  polyhedral  convex  cones 
Cj ,C2,  •  •  •  have  the  sets  •  •  •  ,Tf,,  respectively,  as  (not  necessarily  minimal)  tangential 

representations,  whae  T  =  { « j  ,^2  .•••«,■)• 

The  computer  representation  of  is  generated  by  procedure  INTEJi  using  the  fact  that 

C,+i  =  C,.n{j|x*/i,,.i  2  0)  (3) 

i.e.,  cone  contains  the  vectors  that  are  in  both  cone  Q  and  in  the  halfspace  defined  by  the  plane  perpendicular 


to  vector  and  going  through  the  origin.  Procedure  SOLVE  terminates  as  soon  as  a  cone  of  siu^  POINT  is  found 
in  the  sequence  Cg.Ci  .C2.  ■  ■  ■  .Cj^,  since  all  the  remaining  cones  in  the  sequence  would  also  be  of  shape  point. 

Incorporating  the  shapes  into  the  representations  of  the  cones  simplifies  the  reasoning  needed  to  compute  their 
intersection  with  a  halfspace  because  for  each  shape  of  the  cone,  there  are  only  a  few  possibilities  fcr  the  shape  of 
the  int^section.  Hgure  20  shows  a  state  diagram  in  which  the  nodes  (Le  the  states)  correspond  to  the  shapes  of 
polyhedral  convex  cones  in  R^,  and  the  arcs  (i.e.  the  transitions)  correspond  to  the  possibilities  for  the  shape  of  their 
intersection  with  a  halfspace.  The  iterations  performed  by  procedure  SOLVE  can  be  seen  as  a  sequence  of 
transitions  in  that  state  diagram.  The  initial  state  is  the  shape  SPACE  which  is  the  sh£q)e  of  cone  Cq.  Each  iteration 
in  procedure  SOLVE  causes  a  transition  in  the  state  diagram,  and  the  final  state  is  the  stuqpe  of  the  solution.  The 
actual  transition  is  computed  by  procedure  INTER  which  also  computes  the  necessary  parameters  to  completely 
characterize  the  cone  in  the  syntax  of  figure  18. 

The  inputs  to  procedure  INTER  are  the  computer  representation  of  a  cone  C,  and  a  vector  y  .  The  output  of 
INTER  is  the  computer  rqnesentation  oS  cone  C'  =  C  n  {x  |  x  •  y  ^0).  The  computation  performed  by  INTER 
depends  on  the  shape  of  the  cone  input  Hgure  21  shows  the  cases  in  which  the  shape  of  the  cone  input  is  SPACE  or 
HALFSPACE.  The  Other  cases,  although  more  extensive,  are  not  difficult  to  infer. 

L5.  Example  of  the  Computation  of  the  Directions  of  Feasible  Translations 

Hgure  22  shows  an  assembly  that  has  two  parts  with  seven  planar  contacts  between  them.  The  vectors 
popendicular  to  the  contacts  and  pointing  towards  the  upper  part  are; 

"1=  (1  0  0)  22=  (0  0  1)  23*  (0  2  1)  24=  (0-2  1)  25=  (0  1  1)  26“  (0*1  1)  2?=  (0  0  1) 

For  this  example,  procedure  SOLVE  does  seven  iterations  to  find  out  the  set  of  directions  along  which  the  upper 
part  can  translate.  The  first  iteration  produces  the  intersection  of  the  whole  space  with  the  halfspace  defined  by  the 
plane  perpendicular  to  2 the  intersection  is  the  halfsptux  itself,  whose  representation  is  (HALFSPACE  (1  0  0) ). 

The  second  iteration  produces  the  intersection  of  the  halfspace  obtained  in  the  first  iteration  with  the  halfspace 
defined  by  the  plane  perpendicular  to  n-i,  because  n  j  and  22  are  not  parallel,  the  intersection  has  shape  quadrant, 
and  its  representation  is  (QUADRANT  (1  0  0)  (0  0  1)  ). 

The  third  iteration  produces  the  intersection  of  this  quadrant-shape  cone  with  the  halfspace  defined  by  the  plane 
perpendicular  to  because  2i  >22>  tmd  23  are  linearly  independent,  the  intersection  is  a  polygonal  (triangular) 
cone  whose  rq)resaitation  is  (POLYGONAL  (  (1  0  0)  (0  0  1)  (0  2  1) )(  (0  1 0)  (1  0  0)  (0  -1  2) ) ). 

The  fourth  iteration  produces  the  intersection  of  the  polygonal-shape  cone  obtained  in  the  third  iteration  with  the 
halfspace  defined  by  the  plane  perpendicular  to  24:  because  the  projections  over  24  are  less  than  zero  for  the  first 
edge,  zero  for  the  second  edge,  and  greater  than  zero  for  the  third  edge,  the  rqnesentation  of  the  intersection  is 
(POLYGONAL  ((100)  (0  -2  1)  (0  2  1)  )(  (0  1  2)  (1  0  0)  (0  -1  2)  )  ). 

The  fifth,  the  sixth,  and  the  seventh  itoations  do  not  change  the  polyhedral  convex  cone  produced  in  the  fourth 
iteration,  which  happens  to  lie  entirely  within  the  halfspace  defined  by  the  plane  perpendicular  to  25*  Oie  halfspace 
defined  by  the  plane  perpendicular  to  26>  ^  halfspace  defined  by  the  plane  perpendicular  to  n-j.  This 

conclusion  can  be  made  by  observing  that  the  three  edges  (0  1  2),  (1  0  0),  and  (0-12)  have  greater  than  or  equal  to 
zero  projection  over  25>  26>  ^ 
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Figure  20:  State  diagram  for  procedure  SOLVE 


3: 


procedure  INTER(cone 
case 


(1)  FIRST(cone)  -  SPACE  return  (HALFSPACE  v) 

{2)  FIRST(cone)  =  HALFSPACE  do 

begin 

n  <-  SECOND(cone) 
if  n  and  v  are  parallel  do 

begin 

if  n  and  v  have  the  same  orientation  return  (HALFSPACE  ^ 

return  (PLANE  ^ 

end 


return  (QUADRANT  (v  n)) 
end 


Figure  21:  Part  of  procedure  INTER 


The  firul  solution  returned  by  the  algorithm  is  the  polyhedral  convex  cone  whose  computer  re{»esentation  is 
(POLYGONAL  ((100)  (0  -2  1)  (0  2  1) )( (0  1  2)  (1  0  0)  (0  *1  2)  )  );  this  means  that  any  positive  linear  combination 
of  the  vectors  (0  1  2),  (1  0  0),  and  (0  <1  2)  is  a  feasible  translation  for  the  uppo-  part  in  figure  22.  For  this  example, 
the  result  can  be  verified  by  inspection.  The  set  of  all  directions  d  along  which  translation  is  feasible  can  be  scanned 
systematically  by  letting 

d  =  a  (0  12)  +  b(100)  +  c(0  -12) 

O^a^  1 

O^bsVl-o^ 

C=Vl-fl2-b2 


L6.  Relations  to  Other  Work 

Within  the  research  in  robotic  planning,  the  work  of  Brooks  [7, 8]  has  some  relation  with  the  results  presented  in 
this  paper.  Brooks  formalizes  the  process  of  checking  and  modifying  robot  plans  to  ensure  that  they  will  work  in 
spite  of  inaccuracies  of  mechanical  devices  and  the  inaccuracies  in  the  information  the  robot  has  about  the  position 
and  orioitation  of  parts  within  the  workstation.  That  formalization  leads  to  a  system  of  (not  necessarily  linear) 
inequalities  and  Brordos  uses  a  constraint  manipulation  systnn  to  decide  whether  the  system  has  a  solution  and  to 
find  bounds  for  some  functions  of  the  variables  in  the  plan.  That  constraint  manipulation  system,  however,  does  not 
construct  the  set  of  solutions  to  the  system  of  inequalities.  Moreover,  the  conclusions  drawn  from  that  system  tend 
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Figure  22:  Two  parts  that  have  seven  planar  contacts 


to  be  conservative;  they  ate  safe  to  be  used  in  robot  planning  but  they  may  lead  to  the  elimination  of  plans  that  are 
reliable. 

Systems  of  linear  inequalities  have  been  studied  within  linear  programmming  [36].  where  the  extremal  points  of 
linear  fimcdons  are  sought  in  multidimensional  spaces.  Goldman  and  Tucker  [16]  present  important  theoretical 
results  that  have  been  used  as  basis  for  the  formulation  presented  in  this  paper.  Those  results  alone  have  been  used 
by  Ohwovoriole  and  Roth  [27],  in  the  context  of  mechanical  assembly,  to  solve  a  system  of  inequalities  in  a  five 
dimensional  space.  By  restricting  the  dimension  of  the  space  to  three,  as  we  have  done,  more  efficient  procedures 
could  be  constructed. 

In  addition  to  being  less  efficient  (although  more  general),  the  linear  programming  approach  to  solving  systems  of 
linear  inequalities  has  problems  in  degenerate  cases  which  are  common  in  assembly  planning.  One  degeneracy  is 
the  fact  that  the  set  of  solutions  is  unbounded  in  all  cases,  excq^t  when  no  solution  exists.  Another  degeneracy 
occurs  whenever  the  feasible  solutions  lie  on  a  plane  0-e.  the  set  of  solutions  has  volume  zero);  and  this  happens 
whenevOT  parts  have  parallel  faces  which  are  in  contact  with  other  parts. 
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The  intersection  of  halfspaces  has  been  studied  within  computational  geometry  [30].  Brown  [9]  showed  how  the 
problem  of  finding  the  intersection  of  N  halfspaces  can  be  converted  to  the  problem  of  finding  the  convex  hull  of  N 
points;  that  leads  to  an  algorithm  that  takes  0(N logN)  time. 

Although  more  general,  the  algorithms  in  computational  geometry  have  been  designed  for  the  case  in  which  the 
solution  set  is  bounded.  Like  the  linear  programming  tyjproach,  these  algorithms  have  probl^s  in  degenerate  cases 
which  are  very  common  in  assembly  (banning.  The  procedure  presented  in  this  paper,  has  been  designed  for  the 
assembly  planning  problem.  It  is  less  general  than  those  in  computational  geometry  but  is  more  efficient  since  it 
finds  the  solution  in  at  most  N  steps,  and  it  can  handle  the  degenerate  cases. 

As  mentioned  in  section  1.2,  the  w(^  d  Asada  and  By  [3]  has  some  relation  with  the  results  presented  in  this 
pq)er.  Although  Asada  and  By  modeled  the  contacts  between  the  woikpart  and  the  fixturing  elements  as  point 
contacts,  for  local  translations,  point  contacts  and  planar  contacts  yield  the  same  constraints.  The  ccmditions  that 
Asada  and  By  derive  for  local  translations  are  the  same  conditions  as  equation  2  in  this  pt^.  (The  reado’  is  warned 
that  there  is  an  error  in  equation  20  of  Asada  and  By  pi^)er,  it  should  read  G^-Lq'i.  0.)  But  Asada  and  By  do  not 
address  how  to  determine  the  set  of  solutions  to  equation  2,  which  is  their  equation  20^. 

L7.  Conclusion 

The  problem  of  finding  the  directions  oi  feasible  local  translations  for  a  part  constrained  by  planar  contacts  has 
been  formulated  mathematically  as  that  of  finding  the  set  of  solutums  to  a  system  of  inequalities.  The  system  of 
inequalities  is  represented  by  a  polyhedral  convex  cone,  and  die  solution  procedure  expbits  the  fact  that  in  the  three 
dimensional  space  there  are  only  10  possibilities  for  the  shape  of  a  polyhedral  convex  cone. 

A  syntax  for  the  computer  representadon  of  polyhedral  convex  cones  in  R^,  which  incorporates  explicitly  their 
shapes,  is  presented  along  with  an  implemented  algorithm  that  uses  that  representation  to  produce  the  set  of 
solutirms.  The  algorithm  can  handle  all  possible  cases  and  produces  the  solution  in  at  most  N  (the  number  of  planar 
contacts)  steps.  It  may  take  less  than  N  steps  when  the  only  solution  to  the  system  of  inequalities  is  the  zero  vector. . 

In  addition  to  providing  the  basis  for  testing  the  feasibility  of  assembly  operations,  the  computer  rq>resentation 
generated  by  the  procedure  is  useful  later  in  the  assembly  planning  process  to  guide  the  search  for  a  path,  since  it 
allows  a  systematic  scan  of  all  directions  along  which  local  translation  is  feasible. 


*SiiKe  the  foimulatioa  pretented  by  AsatU  and  By  mclndea  roudona,  the  resulting  system  of  equations  involves  six  variables.  Therefore, 
solving  that  system  would  be  more  complex  than  solving  a  system  of  three  variables  as  Uiat  of  equation  2  in  this  paper. 
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